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Roma,  March  5,  2005 


Workshop  on 

“Construetal  Theory  of  the  Generation  of  Optimal  Flow 

Configurations” 

Roma.  School  of  Engineering,  “Aule  del  Chios  tro”,  17th  and  18th  of  March,  2005 


Constructal  Theory  is  a  new  research  development  aimed  at  the  identification  of  the 
most  efficient  "structure”  that  a  flow  of  mass,  energy,  or  'anything  that  moves7  can  lake. 
The  theory  is  validated  by  comparisons  with  natural  phenomena,  since  it  is  generally 
accepted  that  systems  in  nature  evolve  to  the  most  “efficient"  configuration  under  the 
prescribed  constraints,  Constructal  iheoiy  is  being  considered  for  analysis  of  the 
thermal  management  of  components  and  systems,  for  system  analysis  of  intrinsically 
energy-efficient  vehicles.  For  traffic  management,  for  networks  for  distribution  and 
collection,  etc.  The  approach  shows  potential  for  significant  applications  in  land,  air, 
sea  and  space  systems,  its  generality  being  founded  on  its  theoretical  basis  of  “thinking 
about  design  as  a  science'  The  objective  of  the  workshop  is  to  assess  the  status  and 
limitations  of  present  technologies  for  developing  the  necessary  systems  of  systems, 
current  developments  in  new  technologies;  potential  applications  to  future  system 
developments;  research  directions  for  future  investment. 

The  Workshop  has  been  jointly  organised  by  the  University  of  Roma  1  “La  Sapienza”, 
Duke  University  and  the  United  States  Air  Force  Research  Laboratory  It  has  been 
sponsored  and  funded  by  the  European  Office  of  Aerospace  Research  and  Development, 
Air  Force  Office  of  Scientific  Research  and  by  the  University  of  Roma  1 


The  Organising  Committee 
Adrian  Be] an  (Duke  University) 
David  Moorhouse  (AFOSR) 
Enrico  Sciubba  (U.  of  Romal) 


Program 


Thursday.  March  17  (SnU  Grande  del  Chios tro) 


9.30  - 

Welcome  address 

(R,  Scinbba, 

D,  Moorehouse, 

A.  Bejnn} 

10.00- 

Construct  a!  Theory:  Shape  &  Structure,  from 
engineering  to  nature 

Adrian  iiejan 

1 0.45  - 

Requirements  on  methods  for  system-level  analysis  of 
efficient  Aerospace  flight  vehicles 

David  Moor  house 

1 1.30  - 

Tiee-shiiped  networks  for  distribution,  collection  & 

COOliilft 

Sylvie  Lorente 

12.15- 

Constructs!  Theory  &  variational  principles  for  the 
Navier-Stokes  equations:  is  there  a  link? 

F.nrico  Sciubba 

13.00- 

Lunch  break 

15.00  - 

15.45  - 

Aerosol  particle  deposition,  and  shape  in 
inanimate/ animate  systems 

Antonio  Fr  Miguel 

Flow  architectures  of  lungs  and  river  basins 

A  Heitor  Reis 

16.30- 

Optimal  geometry  for  a  rectangular  slab  with  internal 
convection  cooling 

Adrian  Bejanf 
Mauro  Robbe 

17.15- 

Discussion 

A.Bejan. 

1)  Moorhousc 

|  20.00  -  f  Workshop  Dinner.  Restaurant  “La  Piazzetta* 


Friday.  March  18  (Sala  Piccola  del  Chiostro) 


9.30  - 

Performance  improvement  of  a  tubular  Solid  Oxide  Fuel 
(  ell  .stack  through  the  study  of  its  configurational  design 

Vittorio  Verda 

10.15- 

Numerical  simulations  of  the  Fluid  *&  heat  flow  in  disc- 
shaped  compact  heat  exchangers  with  dendritic  flow 
structure 

Mauro  Robbe 

1  LOO  - 

Energy  usage  in  Aircraft  Equipment  Systems 

Lester  Faleiro 

1 1 .45  - 

On  the  synthesis  of  complex  energy  systems  with  several 
heat  transfer  interactions 

Andrea 

Lazzarctto 

12.30- 

Thermal:  Management  of  Aircraft  Systems 

Dave  Pratt 

13.15- 

Lunch 

15.M1- 

Pane!  discussion 

E.  Sciubba 

15.45- 

General  discussion,  closure  recommendations  for 
future  work 

A.Repn, 

D.  Moor  house 
fvSciubba 

Constructal  Theory: 
Shape  and  Structure, 
from  Engineering  to  Nature 


Adrian  Bejan 

J.  A.  Jones  Distinguished  Professor 
of  Mechanical  Engineering 
Duke  University 
USA 


Shape  and  Structure, 
from  Engineering 
to  Nature 
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Adrian  Bejan  •  Ibrahim  Dincer  •  Sylvie  Lorente 
Antonio  F.  Miguel  •  A.  Heitor  Reis 


1 .  Constructal  theory:  current  progress 


2.  Aircraft  research  and  design:  needs,  current  work 


3 .  Opportunities  for  constructal  theory  in  aircraft  development 


Constructal  theory  ( 1 996) 

Internal  spacings 
Tree-shaped  flow  structures 
Multiple  scales,  nonuniformly  distributed 
Multiple  objectives 

Complex  engineering  flow  systems:  flight 

Transportation,  economics 

Engineering  +  Physics  +  Biology  +  Economics 


Biomimetics 


Nature 


Time 


Engineering 


Constructal  theory 

Nature  — -  Principle 

Time 


Large  resistance  Balanced 

to  fluid  flow  resistances 


Large  resistance 
to  heat  flow 


Optimal  distribution 
of  imperfection 


Optimal 
distribution  of 
imperfection 


Large 

vertical 

resistance 


Balanced 
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Large 

horizontal 

resistance 


Tree-shaped  flow  structure  —  deduced  from  principle 
Round  cross-sections,  ducts  not  seen  in  nature  first 
Cross-sections  of  rivers  not  assumed 

not  ‘simulated’  based  on  an 
assumed  (e.g.,  fractal)  algorithm 


Flow  geometry 


deduced  from  principle 


How,  not  what 


Application 

What 

flows 

How 

Interstices 
high  resistance 
at  the  smallest  (fixed)  scale 

Channels 
low  resistance 
at  larger  scales 

Electronics 

packages 

Heat 

Low-conductivity  substrate 

High-conductivity  inserts 
(blades,  needles) 

River  basins 

Water 

Darcy  flow  through  porous  media 

Rivulets  and  rivers 

Lungs 

Air 

Diffusion  in  alveoli,  tissues 

Bronchial  passages 

Circulatory 

systems 

Blood 

Diffusion  in  capillaries,  tissues 

Blood  vessels,  capillaries,  arteries,  veins 

Turbulent 

flow 

Momentum 

Laminar,  viscous  diffusion 

Streams,  eddies 

Urban  traffic 

People 

Walking  in  urban  structure 

Street  traffic 

Economics 


Goods 


Hand  delivery  and  collection 


Freight,  rail,  truck,  air,  ship 


Constructal  law  (1996): 


“For  a  finite-size  open  system  to  persist  in  time  (to  survive)  it 
must  evolve  in  such  a  way  that  it  provides  easier  and  easier 

access  to  the  currents  that  flow  through  it”, 
or  the  ‘physics’  behind: 


Darwin’s  observations 
fractal  algorithms 


Heat  engine  efficiency 


Constructal  Multi-Scale  Structures 


Ux 


L0 


Package  of  parallel  plates  with  one  spacing. 


Multi-scale  hierarchical  structures 
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T.  Bello-Ochende  and 
A.  Bejan, 

Constructal  multi-scale 
cylinders  in  cross  flow, 


Int.  J.  Heat  and  Mass 
Transfer  48  (2005)  pp. 
1373-1383. 


A.  K.  da  Silva,  S.  Lorente  and  A.  Bejan, 

Optimal  distribution  of  discrete  heat  sources  on  a  plate  with  laminar  forced 
convection, 

Int.  J.  Heat  Mass  Transfer  47  (2004)  2139-2148. 
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A.  K.  da  Silva,  S.  Lorente,  A.  Bejan, 

Constructal  multi-scale  tree-shaped 
heat  exchangers, 


J.  Appl.  Phys.  96(3),  1  August  2004, 
pp.  1709-1718. 
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Minimization  of  work/distance 
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Golden  eagle  • 
Wandering  albatross 
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Ruby-throated 

hummingbird 

►  House  wren 
Bam  swallow 
•  House  sparrow 


Partridge 

•  Pheasant 
Black  vulture 
•  Gannet 

•  Canada  goose 

•  Mute  swan 


Boeing  737 
Airbus  A3 10 
Douglas  DC-10 
Boeing  747 
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Fermat  law 


point-point  flow 


Constructal  law 

volume-point  flow 
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TRADITIONAL  DEVELOPMENT 


SYSTEM  DESIGN  OPTIMIZATION 


Integration  of: 

Airframe  &  propulsion 
+  fully  adaptive  structures 
+  propulsion  &  power  systems 
+  aerodynamics  &  controls 
+  etc,  etc 

Design  for  Synergistic  Interactions 


What  is  the  Common  Measure  for  Comparison  ?? 


Technology  Challenges: 

•  Accurate  prediction/design  tools 

•  Integration  of  airframe  &  propulsion 

•  Control  of  flow  fields  with  plasma 

•  Energy  extraction/power  distribution 

A  Flight  Vehicle  Interacts  Combustion  Energy  +  Mechanical  Work 
+  Thermal  +  Kinetic  +  Potential  Energies  +  +  + 

•1st  Law  Principles  Energy  conservation  only 
•2nd  Law  Principles  ->  Entropy  can  only  be  created 

Minimizing  Exergy  Consumption  Should  Give  an  Efficient  Design 

What  Will  Constructal  Theory  Provide  ????????? 
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Special  Section,  AIAA  Journal  of  Aircraft,  Jan-Feb  2003 


Moorhouse,  “Proposed  System-Level  Multidisciplinary  Analysis  Technique  Based  on  Exergy 
Methods” 

Chandrasekaran  et  al,  “Optimization  of  an  Aircraft  Power  Distribution  Subsystem” 

Paulus  &  Gaggioli,  “Rational  Objective  Functions  for  Vehicles” 

Munoz  &  von  Spakovsky,  “Decomposition  in  Energy  System  Synthesis/Design  Optimization 
for  Stationary  and  Aerospace  Applications” 

Bejan,  “Constructal  Theory:  Tree-Shaped  Flows  and  Energy  Systems  for  Aircraft” 

Figliola  &  Tipton,  “Exergy  Approach  to  Decision-Based  Design  of  Integrated  Thermal  Systems” 

Roth  &  Mavris,  “Generalized  Model  for  Vehicle  Thermodynamic  Loss  Management”  + 

“Method  for  Propulsion  Technology  Impact  Evaluation  via  Thermodynamic  Work  Potential” 


Integration  Framework  for  System-Level  Optimization 
of  Aerospace  Vehicles 

A  flight  vehicle  is  a  ‘device  to  do  work’.  The  complete  set  of  vehicle 
design  requirements  can  be  specified  as  an  energy  system.  The 
mission  can  be  stated  in  terms  of  work  to  be  done  by  the  exergy 
available  from  the  fuel.  The  propulsion  and  every  system  are  then 
components  in  minimizing  the  exergy  consumed.  This  will  provide  the 
necessary  understanding  to  allow  decomposition  into  appropriate 
energy  systems  together  with  appropriate  interactions.  The  exergy- 
based  framework  will  support  the  design  of  a  truly  energy-efficient 
vehicle,  although  a  realistic  flight  vehicle  will  also  be  subject  to  other 
constraints. 


****  * 


•  Define  weight  specific  energy  as  kinetic  +  potential  energies  per  unit  weight: 


E  =  h  +  ^-U2 

2g 


•  Customer  Work  -  includes  generating  the  specific  energy  of  the  payload 
and  overcoming  drag  and  power  requirements: 


dw  dE 

— £-  =  W  — +  PP+DpU 
dt  p  dt  P  P 


REQUIRED 


•  Overhead  Work  -  Sum  of  the  work  consumed  and  drag  caused  by  every 
component  of  the  system  which  has  to  carry  the  payload: 


dw  x— i  (  dE 

— ^  =  y  w, — +d,u 

dt  H  '  dt 


TO  BE  MINIMIZED 


AIAA  Journal  of  Aircraft,  Vol.  40,  No.  1,  Jan/Feb  2003 


The  customer  and  overhead  work  integrated  over 
the  entire  mission  represents  the  the  total  work 
that  has  to  be  supplied  by  the  exergy  of  the  fuel. 

Conservation  of  Energy  in  time  interval  dt  is: 


TTdW 

dw ’ 

dw. 

-tjH  ——  = 

c 

1  o 

dt 

dt 

dt 

H  is  energy  content  of  the  fuel/weight,  r|  is  overall  efficiency,  dW/dt  <  0. 


Design  problem  is  to  minimize  overhead  work  and 
maximize  overall  efficiency 


Design  Mission  Stated  in  Terms  of  Work  to  be  Done 
How  precise  does  this  need  to  be  ? 

Comes  From  the  Exergy  of  the  Fuel  Consumed 

Propulsion  System  Converts  Fuel  Into: 

Mission  Work,  Including  Power  to  Drive  Mission  Equipment 


+ 


Mission  Overhead 

-  Power  to  Lift  Itself  and  Required  Fuel 

-  Overcome  Vehicle  Drag 

-  Power  for  Other  Subsystems 


Waste  due  to  Inefficiencies  in  Operation  &  Thermal  Performance 
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OVERHEAD 

positives 

negatives  * 

Actuation 

essential  functions 

weight  +  power  used 

Structure 

carries  payload 

weight 

Adaptive  structure 

minimize  drag,  provide  control 

weight  +  power  used 

Controls 

trim/  stability/control 

weight  +  cause  drag  +  power  used 

Flow  control 

reduce  drag 

weight  +  power  used  +  waste 

ECS 

essential  cooling 

weight  +  engine  bleed  +  waste 

Propulsion 

supplies  thrust  +  power 

weight  +  waste 

Fuel 

supplies  exergy  +  cooling 

weight 

*  weight  now  reflects  work  that  needs  to  be  done 
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Conventional  high-speed  vehicles  have  a  requirement  to  protect 

the  structure  from  elevated  temperatures  or  provide  active  cooling,  BUT : 


Structures  as  an  Energy  System 

Current  efforts  are  utilizing  structural  deformation  requiring  actuation 
power  as  a  flight  control  effector,  i.e.  wing  twist  for  roll  control. 

In  high-speed  vehicles,  the  structure  will  have  to  be  designed  as  a 
fully-integrated  part  of  the  system.  In  these  vehicles  the  structure 
may  distribute  energy  between  other  systems,  it  may  store  energy 
for  use  as  and  when  required,  and  it  may  be  a  user  of  exergy  for 
tailoring  different  characteristics.  For  the  future,  with  the  advent 
of  active  structures  technology,  there  is  a  need  to  develop  an 
understanding  of  the  vehicle  structure  as  an  energy  system. 


-\dsjp£j vs  SiruciiJjra  Daslg/j 


•  Develop  a  theoretical  framework  to 
identify  energy  flow  inside  of  the  body 
(input  energy,  transferred  energy,  stored 
energy  and  etc.)  for  efficiency  calculation 


Aerodynamic  force 


Actual 
shape  change 


Actuation  force 
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****  * 


Exergy-Based  Design 

Task  required 


In  current  aircraft,  the  engine  obviously  supplies  the  direct  power, 
cooling  air,  etc  for  all  the  aircraft  subsystems.  The  definition  of  a 
turbine  engine  as  an  exergy  component  with  airframe  weight  and  drag 
has  been  done. 

Propulsion  as  an  Integrated  Energy  System 

A  high-speed  vehicle  will  require  new  forms  of  propulsion 
and  airframe  system  integration.  It  may  power  a  device  to  reduce 
airframe  drag  as  opposed  to  providing  thrust  to  overcome  the  drag. 
There  must  be  a  credible  way  to  perform  the  design  trades. 

We  need  to  develop  the  understanding  of  how  the  hypersonic 
propulsion  system  and  all  the  other  required  systems  can  be 
designed  to  common  metrics  and  optimization  criteria. 
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4.  ECS 

X ECS >  Y ECS 
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3.  VC/PAOS 

X  VC /PADS  ’  } T VC  IP  AOS 
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■>  Gy 


VC/PAOS  ■> w VC/PAOS 
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32 


2.  FLS 

X FLS  5  ^FLS 


H 


FLS  ’ KJ  FLS 


u\5 

Weight  of  the  PS 

“35 

VC/PAOS  weight 

“21 

Momentum  drag,  FLS  weight,  power 

“41 

Momentum  drag,  ECS  weight,  bleed  air 

“25 

FLS  weight 

“45 

ECS  weight 

“31 

Momentum  drag,  VC/PAOS  weight,  power 

“51 

Drag,  weight  of  the  AFS 

“32 

Heat  rejection 

“61 

PPAY,  EPAY,  EG  weight 

“34 

Heat  rejection 

“65 

PPAY,  EPAY,  EG  weight 
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A  conventional  flight  control  system  is  a  user  of  exergy  in  the  form  of 
actuation  power,  and  can  impact  vehicle  drag  which  consumes  exergy. 

Future  task  required:  Control  Strategies  for  Energy  Systems 

•  The  control  design  problem  is  magnified  when  we  consider  the 
control  requirements  and  associated  subsystems  of  future  vehicles 

Develop  an  Understanding  of  How  the  Total  Vehicle  Control 
System  can  be  Designed  to  Exergy-Based  Principles. 

Develop  an  Understanding  of  How  to  Control  the  Distribution 
or  Transfer  of  Energy  Between  the  Different  Systems  in  Order  to 
Produce  an  Optimum  Vehicle  Control  System. 


4  +  4  >f  4 


There  is  a  need  to  understand  how  the  results  of  a  laboratory  test 
will  have  to  be  applied  in  the  design  of  a  full-scale  vehicle.  This 
will  become  increasing  critical  as  we  attempt  to  design  a  full-scale 
hypersonic  flight  vehicle  using  components  with  no  previous 
flight  experience. 

Need  to  integrate  mechanical  work,  thermal  energy,  rate  of  change 
of  momentum,  etc,  etc 

This  requires  research  to  develop  an  understanding  of  the 
appropriate  scaling  laws  for  energy  systems. 
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A  common  problem  area  with  flight  vehicle  development  involves 
the  full-scale  characteristics  not  matching  the  predicted  models. 

In  many  areas  of  design  it  is  common  practice  to  design  for 
uncertainties,  e.g.  sensitivity  analyses,  robust  control  theories. 


Uncertainty  Analysis  of  Energy  Systems 

Currently,  there  is  no  database  on  the  expected  accuracy  or  errors 
in  modeling  some  of  the  new  energy  systems  that  will  be  used  in  a 
hypersonic  flight  vehicle.  It  is  probable  that  each  system  will 
interact  with  every  other  one.  The  propagation  of  uncertainties 
through  an  interactive  system  must  also  be  addressed.  There  is  a 
need  to  develop  an  understanding  of  the  uncertainties  in  modeling 
these  systems  and  the  ramifications  of  using  the  models  in  full-scale 
design. 


4  +  4  >f  4 


Optimization  Techniques  for  a  System  of  Energy  Systems 


It  is  not  guaranteed  that  conventional  optimization  techniques  are 
the  most  appropriate  for  the  envisioned  system  of  energy  systems. 
It  is  also  possible  that  a  system  of  energy  systems  will  be  subject  to 
new  and  unforeseen  constraints. 

The  understanding  must  be  developed  of  how  to  minimize  total 
entropy  generation  in  a  complex  integrated  system,  plus  the 
appropriate  constraints. 


Energy  Consideration  is  NOT  New,  BUT  The 
Uses  Have  Been  Implicit,  Disconnected  and  Indirect 
e.g.  the  Breguet  range  equation  ~  work  against  drag 
e.g.  propulsion  system  design 

We  Considered  Explicit  Use  of  Entropy  &  Exergy 
Principles  to  Develop  System-Level  Design/ Analysis 
Methods  Based  on  the  Laws  of  Thermodynamics: 

•  Mission  Requirements  in  Terms  of  Work  to  be  Accomplished 

•  Energy/Exergy  Terms  Flow  Down  to  Every  Component 

•  Computation  of  Entropy  Generation  Rate  ->  Exergy  Destroyed 

•  Minimize  Available  Work  that  is  Wasted,  i.e.  Fuel  Required 


Is  Constructal  Theory  the  Natural  Energy  Flow  for  Mechanical 
Systems  ? 

Does  Natural  Selection  Apply  for  Constrained  Optimization  ?? 
How  do  We  Need  to  Model  Flight  Vehicles  &  Systems  ??? 

What  are  the  Constraints  for  Efficient  Energy  Usage  ???? 

Will  This  Approach  Give  Us  Better  Practical  Solutions  ????? 


What  Further  Research  Is  Required  ? 


T ree-Shaped  Networks  for 
Distribution,  Collection  &  Cooling 

Sylvie  Lorente 

Lab.  Materials  and  Durability  of  Constructions 
National  Institute  of  Applied  Sciences  (INSA), 

Tnu  niiQO  Franro 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations,  Roma. 
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“For  a  finite-size  open  system  to  persist  in 
time  (to  survive)  it  must  evolve  in  such  a 
way  that  it  provides  easier  and  easier 
access  to  the  currents  that  flow  through  it 
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Tree  architectures 

everywhere  in  nature  and  engineering 
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Introduction 


Purpose  - ►  to  make 

a  flow  connection  between 
one  point  and  an  infinity  of  points 


Complex  flow  structures  with  multiple  scales 

:  tree-shaped  networks 

The  flow  resistances  cannot  be 
eliminated.  They  can  be  rearranged, 
assembled  ^ _ 

To  minimize  their  influence  on 
performance  =  optimization 
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method 


integration  of  an  increasing 
number  of  smaller  and  smaller 
flow  components 


elemental  area 


constructs 


Shape  (or  geometry)  is  the  result  of 
the  minimization  of  imperfections,  or 
the  optimal  distribution  of  imperfections. 
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Fractal-like  structure 


T ree  networks  with  loops 
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Networks  at  human  scale 


CTree-shaped  networks  of  insulated  pipes  for  the 
distribution  of  hot  water  uniformly 
over  a  given  territory  (area  A) 


objectives 


Maximum  delivered  temperature 
Minimum  pumping  power 
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Construction  of  the  network 


Elemental  System 


insulating 

material 


r;  inner  radius  of  the  pipe 

r0  outer  radius  of  the 
(pipe  +  insulation) 


Tp,  m0 
AP0 


L0/2 
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Different  configurations 


Case  1 :  every  construct  covers  a  square  area 


L0/2 


AP 


0 
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Case  2:  the  new  construct  is  obtained  by  pairing  the  previous  ones 
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Constraints 


Total  volume  occupied  by  the  ducts 

objective:  minimization  of  the  pressure  drops 

Fluid 

Mechanics  minimization  of  pumping  power 

aspect 


Total  amount  of  insulation  wrapped  around  the  pipes 
objective:  minimization  of  the  heat  losses 


Heat 

Transfer 

aspect 


maximization  of  the  temperature 

of  the  hot  water  received  by  the  end  user 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations,  Roma. 


METHOD 


Pressure  drops  and  heat  losses  are 
minimized  at  every  step 


Set  of  equations: 

Fluid  Mechanics  aspect 


AP  =  C 


L  m2 


r; 


fully  turbulent  flow  (rough) 


T  '  2 

elemental  system  APq  =  q^o_ 5g 


'0 
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example  of  case  1 : 


1  st  construct 


T 


end 


Tv  4m0 

*"a  P., 


AP.,  =  0^° 


(4m0 ): 


+  C'L, 


(2m0  )2  ,  L0  m0 


+  C 


ly\ 


J0 


AP1  =  CL0iho 


12  1 


r 


M 


2r; 


*0 
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total  volume  occupied  by  the  ducts 


0 

-LqTI^  +2L0jir^ 


=  constant 
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Heat  Transfer  aspect 


•  heat  losses  to  the  ambient  (per  unit  of  pipe  length) 


2nk 


In 


V 


vv 


the  pipe  thickness  is  neglected 
k:  thermal  conductivity  of  the  insulating 
material 
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energy  conservation  mCpdT—  qfdx 


T  -  T 

initial  co 


f 

exp 


N  ^ 
InR J 


n  r  _T  2nkL 

with  R  =  -*1  and  N  = - 

mC„ 


"number  of  heat  loss  units" 
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"number  of  heat  loss  units" 


low  N: 

high  performance  materials  (low  k) 
high  mass  flow  rate 


271k  ^ 

2 

•  elemental  system  N0  =  .  n 

m0^p 
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example  of  case  1 


1  st  construct 


T,,  4m 

*~A  P. 


0 


0 


1  - 


T  —  T 

xend  xoo 

Ti-T«, 


=  exp 


N0 

lnR0 


5N0  " 
41nR^ 
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total  amount  of  insulation 


v, 


TiLnr; 


OMq 


3( r, 


v*V 


(r12-1)  +  2(r2-1) 


Fluid  Mechanics 
optimization 


maximization  of  01 


R1  lnR1 
R0  lnR0 


\e) 


1/2 


V  ri1  J 


opt 


Ri  <  Rr 


'opt 


'opt 


The  insulation  shell  is  relatively 
thicker  over  the  smaller  ducts 
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Comparison  of  the  different  systems 


Which  network  serves  the 
farthest  user  better 
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T  <x)  -  =  exp  f-  N  fL  _ ^ _ 1 

T;  -  Tqo  Pl  Jx  xln  (r0  /  rj) J 


Configurations,  Roma. 
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2  levels  of  comparison:  square  areas  4L02  and  16L02 


Tv  4rh0 

A  P1 


"l^i  4m0 

A  P2 


T4  . 

16  m0 

ap4 


the  same 


surface  covered  by  the  network 
amount  of  insulating  material 
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4Ln 


4L?, 


a,,  case  1 


(Tend  ~rroo)"1" 
-p  (Tend  —  Too)  "2" 


amount  of 
insulation 
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a2‘  case  1 


A„,  case  2 


"Case  2"  performs  best 

— ►  Tend  is  higher. 

The  difference  decreases  as  the 
Structure  becomes  more  complex 


insulation 
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The  addition  of  new  users  to  an  existing  network 


Tv  4m0 

A  P., 


lOmn 


_  1  H, 


T 

I 


1 2rh 


« 


one-by-one  »  design 
less  performant  than  the 
constructal  one. 

But,  the  gap  shrinks  for 
larger  structures. 


(a) 


( 


I—  - 


T 

I 


(b) 


14  itiq 


(c) 


1  I  1 
J  _L  L 


I6mn 


(d)  (e) 
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Networks  at  small  scale 

Tree-shaped  networks  in  a  disc-shaped  body 


Micro-channels 
Electronics  cooling 


objectives 


To  deliver  a  fluid  from  a  source 
to  a  given  number  of  outlets 
(users) 

Minimum  AP 


Disc-shaped  area 
Number  of  outlets 
Total  volume  of  the  tubes 
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a 


—  p 
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Flow  resistance  (Hagen-Poiseuille  regime) 

APj  _  128v  Lj 
rh;  71  Df 


When  the  total  volume  of  the  tubes  is  constrained 


Optimal  tube  diameters  follow  Murray’s  law 


Di+1  =  q-1/3 

Di 
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First  step 


AP  =  AP0  +  APj 

L0and  L.,  are  expressed 
using  the  angles  a  and  p 


SAP 

~sp~ 


Leads  to 


p=37.47° 


Finally 


AP 

m 


s  R% 
=  8tw — f 

V2 


With  f  a  dimensionless  resistance 
factor  depending  on  tube  lengths 
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We  obtain  the  connecting  angles 


(48  outlets) 
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Resistance 

factor 
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25 


N  =  constant 

pairing  is  a  useful  feature 

if  N  sufficiently  large 


N  increases 


of  pairing  increases 

:>  complexity  increases 
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25 


By  complexity,  the 
dendritic  flow 
assures  its 
minimal  resistance 

The  designer  can  choose 
between  two  structures 
'  >  ~  same  resistance 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations,  Roma. 


Contents 


Introduction 

Networks  at  human  scales 


Tree-  shaped  networks  of  insulated 
pipes  for  the  distribution  of  hot  water 


Networks  at  small  scales 
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Minimal  pipes  length 


Fractal-like  structure 
Tree  networks  with  loops 
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Comparison:  minimal  length 


•  Every  duct  has  a  fixed  area  allocated  to  it 

•  On  its  area,  the  duct  must  have  the  shortest 
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Minimization  of  flow  path  lengths 


Elemental  curvilinear  rectangle 


Area  of  the  element  fixed 


A  =  ^e(2r8-S2) 


Minimize  L 

L2  =  a2  +  b2 


Result: 


(l-x)2sin 


f 

(2-x)  x  —  1  +  cos 

V 


with  x  =  8/ 


0^ 
2  J 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations,  Roma. 


Construction  of  the  minimal  path  length 


Outer  circle 


Starting  point 

Selection  of  0O:  0.1 


Next 

01  =  20o 

sJ, 

Etc,  etc ... 


Stop  0i  >  2 7t  or  £  <  0 


Note:  rectangle  approximation 
poor  when  r  — ►  0 
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Every  geometric  detail  is  optimized 
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Global  flow  resistance 
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Comparison 


Performance  of  minimal  length  structures 
resemble  closely  the  performance  of  fully 
optimized  structure  Effective  shortcut 
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Every  geometric  detail  is  optimized 

Global  optimization 

The  optimized  design  is  generated 
step  by  step 

Local  optimization 

Assembling  locally  optimized  pieces 
does  not  lead  to  fully  optimized 

structures 


\  X 
.A 
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f =  36.89 

four  pairings,  assumed 
N  =  192 

n0=  12 


Comparison:  fractal-like  structure 


f=  33.16 

three  pairings,  optimized 

N  =  192 

110=24 


f=  23.61 

fom  pairings,  optimized  j  Theory  0f  Generation  of  Optimal  Flow  Configurations,  Roma. 
n0=  12 


Complexity  does  not  mean  performance 
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Tree  networks  with  loops 


Limits  of  the  previous  approach 

In  case  of  accident  along 
the  ducts,  interruption  of 
the  flow. 


Solution:  look  at  nature 
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The  ginkgo  leaf  is  the  exception 
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Designed  networks  mimic  nature: 

The  biggest  networks  have  loops  for  the 
sake  of  security. 


•  Is  an  optimized  network  with  loops  much  less 
efficient  than  a  dendritic  one? 

•  Does  the  result  on  the  comparison  between 
the  two  depend  on  the  complexity  level? 
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Constraints: 

total  surface  (radius  of  the  circle  fixed) 
total  pipes  volumes,  V 


Degrees  of  freedom 

-  number  of  tubes  connected  to  the  center 

-  tubes  diameter  and  lengths 

-  pairing  level 

-  number  of  outlets 
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Laminar  flow 

Local  losses  are  neglected 


If  N  <  5,  optimal  solution:  radial 
network 

N  >  5,  the  «  one  loop  size  »  can 
be  optimized 


very  high  flow  resistances 
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N  >  10,  the  «  two  loop  size  » 
can  be  optimized 


even  higher  flow  resistances 


N-10 

(a) 


N  =  20 
(b) 
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Conclusion: 

The  price  is  heavy  in  terms  of 
fluid  mechanics 
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Other  possibility 


Again,  the  flow  resistances  increase  with  security 
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Robustness  of  optimized 
complex  flow  structures 

Impact  of  having  loops 
decreases  when  complexity 
increases 


Optimal  diameter  ratio  (one  tube  is  cut) 

#  insensitive  to  the  increase  of  complexity 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations,  Roma. 


Concluding  remark 


Several  applications  at  differents  scales 
Urban  hydraulics 
Cooling  of  electronics 


Complexity  is  an  optimization  result. 

Optimized  complexity  must  not  be  confused 
with  maximized  complexity. 
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To  practice  Constructal  Theory  means  “to 
construct  complex  structures  by  connecting 
simple  elements  and  to  optimise  each 
connection ”  [Bejan  2000] 


For  fluid  flow,  “constructing  structures” 
means  to  assemble  flow  paths 

Question:  what  is  the  objective  function  of  the 
optimisation? 
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For  laminar  liquid  flow  in  pipes  &  channels,  constructalists  basically 
contend  that: 

a)  There  is  an  “optimal”  diameter  ratio  at  each  bifurcation,  given  by 
biology-based  “optimal  volume”  considerations  [Murray  1926] 

b)  If  no  heat  transfer  is  involved,  the  “structure”  is  optimised  under  the 
“minimal  flow  resistance”  criterion  (minimal  Aptot) 

c)  If  thermal  transport  is  involved,  the  “maximum  heat  flow”  becomes  a 
suitable  optimising  criterion  (there  may  be  others) 
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My  contention  here  is  that  the  optimising  criterion  is 
always  the  minimum  entropy  generation  compatible 
with  the  prescribed  contraints,  and  that  “Constructal 
Theory”  ought  to  be  viewed  as  an  “entropy-based 
optimisation  tool” 


In  this  first  attempt  to  put  things  in  this  particular 
perspective,  I  will  limit  my  considerations  to  fluid  flow 
only  (no  heat  transfer) 
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1)  The  standard  derivation  of  the  velocity  profile  in  plane 
Poiseuille  flow 

Navier-Stokes  equations 

1  dp 


du  du  du 

- l-u - h  v —  = 

dt  dx  dy  p  dx 

d 


o2  o2  A 

d  u  d  u 


+  V 


V 


dx2  dy 


J 


Physical  assumptions:  steady  flow  (  —  =  0); 

dt  d 


fully  developed  (v=0, 
laminar 

Boundary  conditions:  u=0  at  y=0 

u=0  at  y=d 


=  0  ,  ^P  =  4>  ); 

dx  dx  L 


U  = 


d2  Ap 


, 2 


L  2  ju 


V 


y 

d 


J 
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Once  the  velocity  profile  is  known,  we  can 
compute  the  entropy  generation: 


d2  (Ap)2  ( y2  y|l" 
L2  //T  [d2  d  4y 


And  the  entropy  generation  per  unit  channel 


length 


/IpV  Ld 
L  J  12//T 
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The  pressure  gradient  is  linked  to  the  mass  flow 
rate  by: 

Ap  _  12vm 
L  “  d3 

So  that  the  entropy  generation  per  unit  channel 
length  becomes: 

, .  .  2  T  2t  s,  d  grows  with  L 

_  12vm  L  f  .  m  L  Ld 

*L,d  —  ...  ,5  ,5 

Pl  a  d  sL  d  grows  with  decreasing  d 
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2)  Alternative  derivation  of  the  velocity  profile 
of  plane  Poiseuille  flow 


•  Assume  steady,  laminar,  fully  developed  flow 

*  Assume  that  the  velocity  profile  is  such  that  it  minimises 
the  viscous  entropy  generation  in  the  flow 


s 


V 


unknown  function  of  y 
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The  relevant  variables  in  the  optimisation  are  u,  du/dx=0,  du/dy.  The  Euler- 
Lagrange  equations  for  this  variation  reduce  to: 


d  dsy 


d  dsy 

8y  J 


•  •  • 


g  dsv 

Sy/du^l 

O  — 

UyJ 


The  boundary  conditions  on  u  and  the  additional  constraint  that 
the  flow  must  transport  a  prescribed  mass  flow  rate  provide: 


L  2/^d2  dy 
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This  result  is  indeed  general:  for  viscous  steady  flow  the  Navier- 
Stokes  equations  admit  of  a  variational  formulation  whose 
Lagrangian  is  the  exergy  destruction  (T0*ds/dt).  There  is  no 
definitive  proof  yet  that  this  is  also  the  Lagrangian  for  thermal 
transport  by  convection,  but  it  would  seem  reasonable  to  assume  so. 

The  question  is:  since  the  entropy  generation  “governs”  the 
structure  of  the  flow  in  a  channel,  why  should  it  not  be  the 
“Lagrangian”  for  the  structure  of  more  complex  flow  domains? 


I  shall  now  try  to  derive  an  “optimal”  bifurcation  structure  from  a 
“ minimum  ds/dt ”  principle 
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The  relevant  variables  are:  m0,  mv  L0,  L1?  d0, 
dj,  and  the  fluid  properties. 


a)  m0  is  supposed  prescribed 

b)  nij  =  0.5  m0 

c)  The  flow  is  isothermal 


EoU  2 

do 


miLi 


moLo 


1  + 


Lidp 

LodfJ 


We  want  to  minimise  stftt  w.r.t.  Lft,  Lt  and  to 
the  ratio  dj/dQ 


li 
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In  general,  for  a  prescribed  global  geometry,  L0  and  Lj  will  be  functionally 
linked: 

Lrf(Lo> 

It  follows: 


2t  f 


A  _  k  m0L0 

»tot  -  ^ 


f  d  l  V 


\ 


L0  + 


f(L0) 

S5  ) 


where  8=d1/d0 


which  has  non-trivial  minima  that  depend  on  f: 


ds 


tot 


dL 


=  K 


m2 
1,1 0 


0 


'  d^  V 


1  +  ^1 


S5  ) 


b  +  S* 


'0,opt 


and  LQ?0p^  — 


if  f=aL02+bL0+c; 


a 


ink!) 


aln(a) 


if f=exp(-aL0) 
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If  the  link  between  L0  and  Lx  is  one  of  pure  proportionality,  Lj=aL0: 


dstot  =Kmlf 


dL 


0 


f  ~j5 


1  + 


a 


\ 


d;JV  SD) 


and  no  minimum  is  possible 


But  if  the  link  between  L0  and  is  dictated  by  a  conservation  of  the 
total  volume  occupied  by  the  channels, 

I  _  Vtot  L0d0 
1  2d!  2di 

it  turns  out  that  there  is  a  minimum  for 


v 


d 


0  ^opt 


1 

w 
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In  conclusion,  there  is  a  set  of  different  minima,  each  one 
corresponding  to  a  proper  (i.e.,  mathematically  consistant 
and  physically  realistic)  link  between  L0  and  L, 


Let  us  consider  an  application  to  a  small  set  of 
bifurcational  geometries:  which  one  is  “optimal”? 
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5  => 

A’  1 
0.891  J 

i  A' 

L  0.794 

B, 

0.891 

B, 

0.794 

c'  / 
0.891  / 

V 

\0.794 

d  j 

0.891/1 

D, 

\0.794 

o 

H-l 

Si/ 

/4 

Si 

0 

uS 

Li 

1/ 

/i 

i 

uS 

c 

tot,adim 


1.483 


tot,adim 


1.611 


Ap 


tot,adim 


1.483 


1.09 


Time 

through 

L0+L, 


1.611 


1.483 
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These  results  -albeit  grossly  preliminary-  indicate  that: 

1)  “minimum  Ap”  and  “minimum  s”  generate  the  same  optima 

2)  “minimum  volume”  consistently  indicates  different  optima 

3)  “minimum  delivery  time”  leads  to  the  same  optima  as  the 
minimum  volume 

(1)  was  to  be  expected:  under  the  present  model,  Ap  oc  stot 

(2)  the  reason  for  which  minimum  V  and  minimum  t  lead  to  the  same 
optimum,  remains  to  be  investigated 

(3)  biological  scientists  generally  assume  that  “minimum  V”  is  the 
objective  function  for  the  optimisation  of  dendritic  structures,  but 
it  is  more  likely  that  Nature  may  prefer  something  like: 
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£=  aL0,d0Stot  +  /?L0,d0Vtot 


I  II 

i.e.,  an  “exergo-economic  criterion! 

(Suggestion)  “Natural”  structures  are  generated  under  the 
following  conditions: 

1)  A  weighted  sum  of  the  exergy  destruction  (friction,  thermal 

gradients...)  and  of  the  embodied  exergy  (material  &  energy 
to  “build”  the  channels)  is  minimised; 

2)  External  boundary  conditions  dictate  the  details 
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L=  ffL0,d-s 


o 


tot  P  L  o  o  ^  tot 


II 


For  instance,  for  the  “constructor”  of  the  bifurcations  discussed  above: 


1)  If  the  pumping  power  required  to  make  up  for  the  entropy  dissipation  is 
much  smaller  than  the  exergy  needed  to  build  the  channels,  a«p,  and 
the  optimal  structure  is  Cl; 

2)  If  exergy  is  abundant,  and  p«a,  the  optimal  structure  is  A2; 

3)  If  a=P=0.5,  the  optimal  structure  is  D1 
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L=  a  L  o  ,d  o  S 


tot  P  L  o  o  ^  tot 


I 


Constructal  Theory  &  Exergy  Analysis  have  the  same  goal: 
are  they  using  the  same  tools? 

There  is  an  enormous  incentive  to  investigate  this  field,  for 
the  applications  in  engineering,  physical  &  biological 
sciences  are  of  substantial  importance  (heat  transfer,  filters, 
chemical  reactors,  composite  structures,  allocation  of 
resources,  neuronal  science,  vascular  systems...) 
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Conclusion: 


Constructal  Theory 


is  good  for  the  scientist : 
wnetner  it  is  good  for  the  rest  of  mankind , 
it  is  arguable ” 


Erwin  Chargaff,  1905-2002 
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A  set  of  inanimate  particles 

coagulation  =  agglomeration 


Deposition 

3 


Katzer  et  al.  2001 


8 


Smith  et  al.  2002 
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A  set  of  living  organisms 


(stinging  cells)  i 


zooxanthellae 


living  tissue 
linking  polyps 


nematocyst 


Anatomy  of  a  Coral  Polyp 


wmm 

Wwm&m 


tentacles  with 


nematocysts 


r 


Coral  colony  association  of  polyps 

Tentacles  to  trap  food 

Cells  on  the  lower  sides  and  bottom  produce  the 
limestone  that  builds  the  colony 


/  f  / 

gastnovascularcavity  / 

\  (digestive  sac)  skeleton 

V  ..  '  . 


limestone  calice 


Kaandorp  and  Sloot 
(2001) 
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Why  agglomeration/deposition  is  important  ? 


•  medical  care  for  people  with  chronic  respiratory 
diseases  cost  45.3  billion  USD 


European  Respiratory  Journal  (2002) 


inhalation  thorapy  (drug  delivery  via  the  respiratory  system  easier 
than  delivery  via  ingestion  and  the  transdermal  route) 


Air  pollution  control  devices 


•  Electrostatic  precipitators  (expansive,  fail  to  filtrate  fine  particles  and 

high-resistivity  duSt)  Dietz  (1977),  Shapiro  (1998) 

•  Space-charged  precipitators  (only  effective  for  high  density  aerosols) 

Shapiro  &  Brenner  (1990) 

•  Cyclones  (only  effective  for  coarse  particles  >  5  pm)  Shapiro  &  Brenner  (1990) 

•  Diffusion  batteries  (very  efficient  for  particles  <0.1  pm)  Knutson  (1999) 

•  Filters  (simple,  highly  efficient  for  both  micrometer  and  submicrometer 

particles,  coast-saving  process)  Shapiro  &  Brenner  (1990),  Miguel  (2003) 
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Outline 

Constructal  view 

Particle  coagulation  (agglomeration) 
Colony  association  of  polyps 
•  Design  of  air  pollution  control  devices 
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Design  of  air  pollution  control  devices 


C)  Operate  at  very  low  Re  (laminar  flow) 

Very  important  for  removal  of  submicrometer  particles 


What  are  the  rules  to  obtain  the  best  design  ? 

Constructal  idea: 

The  ultimate  purpose  is  to  bring  the  system  into  equilibrium 
trough  capturing  the  all  particles  in  the  air. 


Global  constrain 

Due  to  practical  and  economical  reasons  the  device  should  fit  into  a  volume  V 


i=> 

c  ^ 

''■'in 

U 


■=> 


\r 


4>=H7d 


Similar  to  optimal  geometry  of  plates  in  forced  convection  heat  transfer 

(i)  finding  of  two  d-value  limits  Bejan  <20#0);  Bejan  <20#4) 

(ii)  identifying  the  d-value  that  maximizes  the  particle  deposition 
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(i)  d  (small)  ^  (large) 


K 


-H 


c  ^ 

U 


-d 


H 


outlet  concentration  of  particles  (minimum) 

steady-state  particle  transfer  density  m  -  (q 

p  HLW  1 

fully  developed  (Hagen-Poiseuille)  flowy  _  d  Ap 

12p  L 


-C0) 
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d  (large)  <|>  (samll) 


steady-state  particle  transfer  density  for  the  boundary  layers 


2JXLW 
p  HLW 


m  =  L238^^— 

HL  (  p  ; 


(cm-c0) 


Sh  =  0.664  Re"2  Sc1'3 

Sherwood  number  (AIVD) 
Schmidt  number  (o/D) 


free-stream  velocity  U0 

ill .5  _  ApH 

°  1 .328Ll'5(pp)0'5 


Bejan  (2000,  2004) 


-Nu 
— Pr 
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(ii)  optimal  number  of  plates  method  of  intersecting  asymptotes 


opt 


L 


=  2.73Be 


-1/4 


m 


Bem  = 


L2Ap 


Bhattacharjee  &  Grosshandler  (1988) 

Best  elementary  construct 


K 


fully  developed 
region 
K - 


flow  fluid  developed  particles  interact  with  plate 
air  in  the  core  duct  not  interact  with  plate 
Maximum  particle  transfer  L^ev=L0pj 


(  d^l 

=  1.86Sh1/4Bem1/4 

IlJ 

opt 

-!Jdev<IJ 
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Device  with  a  porous  structure 

(i)  finding  of  two  extremes 

£  ~  0 


outlet  concentration  of  particles  (minimum) 


steady-state  particle  transfer  density  mp 


U 

L 


(cin-cj 


Darcy  flow  u  =  — —  Carman-Kozeny  equation  k  =  c  -  d2 


p  L 


(l-8f 
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8  ~  1 


steady-state  particle  transfer  density  to  N  collectors 

m  =  Nc  D  — 

p  T>m  ^ 


(Cin-CD) 
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(ii)  identifying  the  S-value  that  maximizes  the  particle  deposition 


particles  transferred 
from  air  stream  to  V 


particles  transferred 

from  air  stream  to  collectors  of  porous  material 


<;K 


p  (1  -  e)2  L2 


■Ap(Cm-C0) 


?ms(l-s)|-(Cm-C0) 
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Constructal  view  of  aerosol  agglomerate  on/deposM on  process 


What  determines  the  shape  of  in  these  processes  ? 


Constructal  law  .  ►  architecture  optimized 

'  '*■  accumulation  of  particles  maximized 


aerosol  particles  >100  nm  carry  an  electrostatic  charge 

Vincent  (1986);  Kerasev  et  al.  (2001) 


aerosols  travel  under  the  influence  of  Coulomb  forces 


t=0  t>0 

Electric  field  from  a  point  charge  Coulomb  Force  AerosolO  migrates  Stokes  flow  velocity  U=dr/dt 

E=kq/r2  F=qE=kqq/r2  0  |_id . . 

F  =  3ti  £ — U 
Cc 

tx  =  — x3  *  “flying”  time 

Cck  pd 


The  flow  of  particles  that  lands  on  0  originates 

a  spherical  annulus  of  radius  x  and  thickness  Ax  AN  =  n  x2  Ax 


AN/At  ~  d(D3)/dt 


♦  D  = 


2kccn 


•  Cluster  diameter  grows  fast  early,  and  slow  late 
•Mass  accumulated  ^D3)  increases  linear  with  t 


pd 

qq 


xl/3 

t 

J 
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C)  Constructal  question 

r* 

There  is  any  another  configuration  (shape)  that  maximizes  the  mass  accumulated  i 


The  oldest  site  The  newest  site 


Needle  body  of  revolution  with  radius  R=fx1/3 


x  instantaneous  needle  length  (Ut) 


2kccnpd 

uqq" 


\  1  /  3 


V 


J 
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V, 


sphere ~ 


4  Da 
2 


v^-y 


4 

3 


X 


—  — Tcf^X  ^needle  — 

0 


3tt  ,2  5/3 


x=Ut 


Aggregate 

volume 


Largest  sphere 
before  needles  — ► 


0 
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zooxanthellae 


rr tpy 


living  tissue 
linking  polyps 


nematocyst 


Anatomy  of  a  Coral  Polyp 


tentacles  with 
nematccysts 
(stinging  cells)  ^ 


mouth 


/  v-  * 

/  w 

gastrovascular  cavity 
\  (digestive  sac) 

v.  \ 

limestone  calice 


Constructal  view  of  coral  colony 
association  of  polyps 


Tentacles  to  trap  food 


*  Cells  on  the  lower  sides  and  bottom  produce  the 
limestone  that  builds  the  colony 


skeleton 


©  1999  Encyclopaedia  Britannica,  Inc. 


Porites  spp.  (stony  coral) 


rounded  and  compact  dendritic 


cm/year 


cm/year 


growth  rate  1.2+0. 5  0.2  to  4.8 


Lough  and  Bames  (1997) 


Harriott  (1999) 


Constructal  view  of  coral  growing  process 

Why  rounded  and  compact  forms?  Why  branches? 
Why  are  the  branches  necessary? 

Food  particles  very  small  convective-dispersive  transport 
Growth  process  starts  at  birth  (t=0)  T  _  t  t  t 

-^C  ^  C  1 

Sheltered  growth  site  /  „  . 


Pe  = —  <1 
l  D  *) 


food  particles  near  colony  Cc  (<cbulk) 


food  “wave” 


L  p-(D*t) 


Food  particles  velocity 


fast  early  and  slow  late 
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How  can  the  access  to  nutrients  by  the  colony  of  coral  be  optimized  ? 


U„  ~  const 

v 


fast  early  and  slow  late 
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Does  the  nutrient  size  influence  the  form  of  coral  colony  ? 


Cc„kBT 


D* 

71 


U„  ~  const 


Colony  more  exposed  to  water  current 


Pe 


un 


D* 


>1 


) 


Uc  1 .2+0.5  cm/year 


Lough  and  Bames  (1997) 


0.2  to  4.8  cm/year  Harriott  (1999) 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations 

Universita  Degli  Studi  di  Roma 


What  is  the  constructal  view  of  stony  coral  growth? 


Branches  generation  allows  a  faster  growth  of  the  coral  colony 
because  it  provides  internal  paths  such  that  the  entire  substrate 
volume  approaches  internal  equilibrium  faster  (depletion  of 

nutrients). 


Bacteria  colony  &  plant  roots 


Miguel  (2003) 
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Final  remarks 


Constructal  Theory 


•  design  optimization  of  air-cleaning  devices 

•  cluster  of  particles  grows  fast  early,  and  slow  late 

•  at  small  times  spherical  agglomeration;  at  long  times 
needle  (dendrite)  more  effective 

•  in  stony  coral  branches  generation  provides  the  best 
internal  path  to  reach  the  food  (i.e.,  internal  paths  such  that 
the  entire  substrate  volume  approaches  internal  equilibrium 
faster). 


Workshop  on  Constructal  Theory  of  the  Generation  of  Optimal  Flow  Configurations 

Universita  Degli  Studi  di  Roma 


Workshop  J-'C  onshuc  tal  Theoty  of  (he  Generation  of  Optimal  How  C  onfigwo  Hon”,  Roma ,  17/19  Match  200 5 


Construe rill  iiisoxy  o£  flow 

iii'cliiiociiii'Scj  oi  iii£  luilfifS 


(Reis,  Miguel  and  Ay  din,  2004) 
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Outline 


•  A  fluid  i/33  with  purpotd 

•  Bronchia]  ^aa  resistance  and  ali/aolarresisflanc 


j  Optir/ktion  okihe  respire  toiytraa  bs j 
ConSl/J.J£tEjJ  RflftCiple 

j  ConstaictaJ  lawisofihe  Ju/j^j  ainNays 

J  C  £)/]£  JiJ£J ajJ3 


d  on  lIj 


h  2C05 
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The  respiratory  tree 

•  Starts  at  the  trachea; 

•  Channels  bifurcate  23  times 
before  reaching  the  alveolar  sac. 


Has  this  special  flow  architecture  been  developed  by 
chance  or  does  it  represent  the  optimum  structure  for  the 
lung’s  purpose,  which  is  the  oxygenation  of  the  blood? 
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The  purpose  of  the  respiratory  tree 

•  Oxygenation  of  blood; 

•  Removal  of  carbon  dioxide. 


The  trade-off 

•  High  alveolar  surface  promotes  better  oxygenation,  but  requires 
increased  access  to  the  external  air  (higher  flow  resistance). 

•  Lower  level  of  channel  branching  reduces  flow  resistance,  but 
reduces  also  the  number  of  alveoli  and  therefore  decreases  the 
area  allocated  for  blood  oxygenation. 
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The  asymptotic  leeways 

I)  A  duct  system  that  ends  with  an  alveolar  volume  from  which  the 
oxygen  diffuses  onto  the  tissues,  where  it  meets  the  blood,  and 
in  which  the  carbon  dioxide  diffuses  after  being  released  from 
the  blood; 

tFl0w=  rrr|L2/(d2AP)~l  s 

II)  A  unique  volume  open  to  the  external  air,  in  which  the  oxygen 
reaches  the  blood  in  the  tissues,  and  removes  the  carbon 
dioxide  rejected  from  the  blood,  only  by  diffusion  through  the 
internal  air. 

t[Mf  =L2/D-104  s 
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Bronchia!  tree  resistance 

Assumptions 

•  Airflow  within  the  bronchial  tree  is  assumed  to  be  laminar 
(Re<1000),  isothermal  and  incompressible 

•  The  difference  in  chemical  potential  drives  oxygen  into  the  lungs 
and  cabon  dioxide  out  of  the  lungs  Ap  =  p-1AP  +  As 

•  The  bronchial  tree  is  composed  of  cylindrical  channels  with 
Hagen-Poiseuille  flow, 
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Airflow  rate: 
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Minimum  flow  resistance  at  a  bifurcation 


Constructal  rules : 


D  ID  =2'1/3 

unun-l 


L  II  =2'1/3 

hrS-l 


diffusive 

tesistance 


Resistance  posed  by  the  nth  bronchial  tube 
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^bbn  _  _  ^sbn 
rbn  rbn 


Resistance  due  to  bifurcation 


rbn 


m0  on/3 
87ip2Do 


Resistance  of  the  nth  level 
(2n  bronchial  tubes) 


128vLo  (,  mo  2_2n  1 3  A 
rcpDj  t  +  1024pvL0  , 


overall  convective 
resistance 


„  K  128vLn 

P  — \  r  —  0 

rH+m0(l-2-™»)l 

b  /  j  n  -p.4 
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m0x  =  ^  )o  “  tax  hb  = 

^  VKox^b 


<|>  -  concentration 


mox=2 


N 


2*dPDox(AlOa 


RoxT 


N 

d  =  L-£Li 

i=l 


d  =  4.85x2“(N+1)/3L 


0 


resistance  of  the  Nth  level  of  bifurcation 
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D  rv 

^OX  — 


256vLq 


((^OX  )()  ^OX  ^P 


(  N +1  } 


0.13x2 


-2N/3 


K)  T 

v  s  7OX 


ttLqD^oxP 


Air,  oxygen  and  carbon  dioxide  properties  were  taken  at  36°  C 
<(>ox  is  evaluated  from  the  alveolar  air  equation  :  ((4>ox)0_<t>ox)Q_S=0, 

Q  is  the  tidal  airflow  and  S  is  the  rate  of  oxygen  consumption. 
((()ox)air=0.2095;  (c|>cd)air~0.315xl0‘3;  Q~6  dm3/min  and  S-0.3  dm3/min. 

trachea  length  L0~  15  cm;  trachea  diameter,  D0  ~1.5  cm' 
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Constructal  law  (principle)  the  flow  architectures  evolve  in  time  in  order 
to  maximize  the  flow  access  under  the  constraints  posed  to  the  flow 


2. 164  In 
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OX 
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Minimization  of  the  global  resistance  to  oxygen  access  with  respect  to 

level  of  bifurcation 


Nopt  =  23  (23.4) 

Minimization  of  the  global  resistance  to  carbon  dioxide  removal  with 

respect  to  level  of  bifurcation 


Nopt  =  23  (23.2) 


Resistance  (Js/kg2) 
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Total  resistance  to  oxygen  and  carbon  dioxide  transport  between 
the  entrance  of  the  trachea  and  the  alveolar  surface  is  plotted  as 

function  of  the  level  of  bifurcation  (N). 


Level  of  bifurcation 
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Constructal  laws  of  the  human  respiratory  tree 


(I) 

If  the  number  N  =23  is  common  to  mankind  then  a  constructal 
rule  emerges: 


“the  ratio  between  the  square  of  the  trachea  diameter  and  its 
length  is  constant  and  a  length  characteristic  of  mankind:" 

D2 

X= — -=const  .=1 .6x1 03  m 
Lo 
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V  -  Global  volume  of  the  respiratory  tree;  L  -  average  length 
A  -  Area  allocated  to  gas  (02  and  C02)  exchange 

(  _  .  A 1/2 


— =8.63AL 
Lo 


V 


vD0X(|)0X 


L(R«  L  Tfe+ox )«  >) J 


oil 

The  non-dimensional  number  ALA/,  determines  the  characteristic  length 
X=D02ll0,  which  determines  the  number  of  bifurcations  of  the  respiratory  tree 
by: 

“The  alveolar  area  required  for  gas  exchange,  A,  the  volume  allocated  to 
the  respiratory  system,  V,  and  the  length  of  the  respiratory  tree,  L,  which 
are  constraints  posed  to  the  respiratory  process  determine  univocally  the 
structure  of  the  lungs,  namely  the  bifurcation  level  of  the  bronchial  tree.  ” 
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Global  volume  and  average  length  of  the  respiratory  tree,  and 
area  allocated  to  gas  (02  and  C02)  exchange 


Average  length  of  the  respiratory  tree:  L=4.85L0=0. 68m  ; 

Observed  ~lm 

Alveolar  diameter:  d=6V/A=2.67xlO'3m  ;  Observed  ~  10  3m 
Alveolar  surface  area  as  A=2237id2=187.7m2 ;  Observed  100  -150m2 


Total  alveolar  volume  V=223(7i/6)d3=83.5dm3 ;  Observed  8  -10m3 
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CONCLUSIONS 


•  The  best  oxygen  access  to  the  tissues  where  it  reaches  the  blood  is 
performed  by  a  flow  structure  composed  of  ducts  with  23  levels  of 
bifurcation 

•  The  same  structure  has  been  shown  being  optimized  for  carbon 
dioxide  removal. 

•  The  optimized  number  of  bifurcation  levels  matches  the  23  levels  that 
the  physiology  literature  indicates  for  the  bronchial  tree. 

•  Theory  also  predicts  the  dimension  of  the  alveolar  sac,  the  total 
alveolar  surface  area,  the  total  alveolar  volume,  and  the  total  length  of 
the  airways.  These  values  agree,  at  least  in  an  order  of  magnitude  sense 
with  the  values  found  in  the  physiology  literature 
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•  It  was  shown  that  the  length  A  (defined  as  the  ratio  between  the  square 
of  the  first  airway  diameter  and  its  length)  is  constant  for  every  individual 
of  the  same  species  and  related  to  the  characteristics  of  the  space 
allocated  for  the  respiratory  process. 

•  The  length  A  is  univocally  determined  by  a  non-dimensional  number, 
AL/V,  which  involves  the  characteristics  of  the  space  allocated  to  the 
respiratory  system,  namely  the  total  alveolar  area,  A,  the  total  volume  V, 
and  the  total  length  of  the  airways,  L. 

•  In  general,  we  conclude,  that  for  every  species  whose  respiratory  tree  is 
optimized  the  same  rule  must  hold,  and  exhibit  the  respective 
characteristic  length,  A. 


•The  application  of  the  Constructal  Principle  to  the  generation  of  the 
optimal  configuration  of  the  respiratory  tree  was  based  on  the  view  that 
Nature  has  optimized  the  living  flow  structures,  in  time. 
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Naiiue]  flow  architec  ts 
£2  sj  Jj  jj£j  J;j  ws  of  it t/er  ba sins 
River  networks  33  c  onsfxyctaJ  fluid  1/333 
A  constructs]  xxiodel  of  river  basin  development 
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Tagus  river  basin  (Tejo,  in  Portuguese) 


Length  (km):  1  100 
Drainage  Area  (km2):  88  700 
Discharge  (m3/s):  600 
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Tevere  river  basin 

Length  (km):  405 
Drainage  Area  (km2):  17375 
Discharge  (m3/s):  250 
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River  networks  are  self-similar  structures  over  a  range  of  scales 


Horton’s  law  of  stream  lengths:  Lf  / LMj  - 


Hack’s  law  La=  a(AJb 

b  -0.56  -0.5 


Melton’s  law:  Fs  =  0. 694(  L)0) )~ 

D„  =  L/A  ■  F.=  NJA  s 


1.5  <Rl<3.5 


Horton’s  law  of  stream  numbers:  ni-llni  =RB 
3  <Rb<5 


00-3 
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CONSTRUCTAL  THEORY  (Bejan) 
River  basins  as  area-to-point  flows 
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First  construct  made  of  elemental  areas,  A0  =  HJL0. 

A  new  channel  of  higher  permeability  collects  flow  from  the  elemental  areas. 
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The  optimised  geometry  of  area-to-point  flow 
(channels  with  Hagen-Poiseuille  flow,  Bejan) 

K  =  K/A0  (Hi,Li)  =  (Hi,Li)/(A0)172  &i=Di/Hi 
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Void-allocation  (channel)  optimisation 
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Assumption: 

Channel  hierarchy  is 
understood  in  the  Hortonian 
sense,  i.e.  all  streams  of 
order  i  are  tributaries  of 
streams  of  order  i+1. 
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Horton’s  law  of  stream  lengths 
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Number  of  streams  of  order  i :  Nj  =  Jif  x  Tlj_j  xrif_2 


The  ratio  of  the  number  of  streams 
of  order  l-l  to  the  number  of  streams 
of  order  /  is  Nf  / Nf+ 7  =  n\ 
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Horton’s  law  of  stream  numbers 


Nf/Nj+ 7  -{^>qK  ^  ^  j  ~  constant 
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Relationship  between  RL  and  RB 


Horton’s  law  of  stream  lengths: 
1.5  <Rl<3.5  ;  Rl  ~2 

Li!Li-l  =  RL 


Horton’s  law  of  stream  numbers: 

3  <Rb<5  ;  Rb  ~4 


ni-l!ni=RB 


-1/3 


*fi/Ni+I 


Rb  ~  (Rl)2 
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Hack’s  law 

K  =  a(Af 

b -0.56 -0.5 

A,  =  HtL, 


Hi 

Lt 

ni 

0 

25/63l/6£l/6 

&y2 

~ 

0,1/2 

9 0 

25/63l/6£l/6 

1 

21/6&‘/2 

&j/2 

<P3/2<P*/2 

3!  /  6  £1/6 

23  /  2  £1/  2 

24  /  3 3I/  6  £2  /  3 

2 

&j/2 

31/6(020,)3/2 

3i/3(020  )322 

2i/2K,/2 

25/301o/2K5/6 

25/6<P0K223 

3 

31 26  (<P2&,  )3/  2 

31 2  3  (<P3<P2&,  )3  2  2 

21  /  6  3I  /  3  (0302f  2  2 

22/34>I0/2k5/6 

24/30oK7/6 

<P0K2  2  3 

4 

32/3(&3&2®i)3/2 

3,23(0403020j)322 

27  26 322 3 (®4V>3)3  2 2 

21/30oK7/6 

2,/203o22K9/6 

0ok223 

Constructal  Hack’s  exponent  /?for  river  basins  up  to  order  4. 


, 0.750 
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constructal  rule  for  the  exponent  ib 


+ 1 
4co 


Muller  reported  that  b~  0.6  for 
river  basins  less  than  8000  mi2,  b 
~ 0.5  for  basins  between  8000  and 
100,000  mi2  and  b  ~  0.47  for 
basins  larger  than  100,000  mi2 
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Melton’s  law:  Fs  —  0.694{DCO)^ 

D„  =  LJA  ;  F.=  NJA 
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A  constructal  model  of  river  basin  development 

(Errera  and  Bejan) 

Objective:  providing  minimal  resistance  to  flow 

Changes  in  the  river  channel  are  possible  because  finite  blocks  can 
be  dislodged  and  entrained  in  the  stream 

Rule:  every  squared  block  of  area  L2  and  height  W,  is  dislodged 
whenever  the  pressure  difference  AP  across  the  block  surpasses  the 
critical  force  needed  to  dislodge  it,  i.e.  APLW  >  xL2, 

Flow  resistance  decreases  with  block  removal.  Further  increase  in 
the  flow  rate  may  create  the  conditions  for  the  removal  a  second 
block  and  for  the  repetition  of  the  process.  A  macroscopic  dendrite¬ 
like  structure  emerges  progressively  while  flow  resistance 
decreases. 
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Conclusions 


The  scaling  laws  of  geometric  features  of  river  basins  can  be  anticipated 
based  on  Constructal  Theory,  which  views  the  pathways  by  which 
drainage  networks  develop  in  a  basin  not  as  the  result  of  chance  but  as 
flow  architectures  that  originate  naturally  as  the  result  of  minimization  of 
the  overall  resistance  to  flow  (Constructal  Law). 


The  ratios  of  constructal  lengths  of  consecutive  streams  match  Horton’s 

law  for  the  same  ratio. 


The  ratios  of  constructal  lengths  of  consecutive  streams  match  Horton’s 

law  for  the  same  ratio. 


Agreement  is  also  found  with  the  number  of  consecutive  streams  that 
match  Horton’s  law  of  ratios  of  consecutive  stream  numbers 

Hack’s  law  is  also  correctly  anticipated  by  Constructal  Theory,  which 

provides  Hack’s  exponent  accurately. 

Melton’s  law  is  verified  approximately  by  Constructal  Theory,  which 
indicates  2.45  instead  of  2  for  the  Melton’s  exponent. 

It  was  also  shown  that  a  constructal  model  of  erosion  can  generate 
dendrite  like  patterns  of  low  resistance  channels  that  reproduce  the 

typical  shape  of  a  river  basin 

As  has  been  demonstrated  with  many  other  examples  either  from 
engineering  or  from  animate  structures,  the  development  of  flow 
architectures  is  governed  by  the  Constructal  Law. 


Optimal  geometry  for  a 
rectangular  slab  with  internal 
convection  cooling 


Initial  Assumptions 


The  slab  thickness  is  prescribed 

The  geometry  is  made  dimensionless  w.r.t.  the  slab 
thickness 

Some  characteristic  length  ratios  between  relevant 
structures  are  also  prescribed 

Additional  geometrical  constraints  are  introduced  to  make 
the  model  as  real  as  possible  and  as  close  as  possible  to  the 
sections  of  internally  cooled  blades 

The  boundary  conditions  are  realistic  and  derived  from  real 
blade  cooling  cases 


Slab  configuration 


Used  relations: 


—  *D  2 

4  0 


Corrected  relations: 
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Values  of  surface  heat  transfer  coefficients 
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STATIC  TEMPERATURE 
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Exch.Coeff. 


Dimensions 


fil  value 

fi2  value 

fi3  value 
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Diameter  H2 
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|  Contours  of  Static  Temperature  (k) 

Mar  14, 2005  Contours  of  Static  Temperature  (k) 

FLUENT  6.0  (2d,  coupled  imp,  lam) 

Mar  14, 2005 
FLUENT  6.0  (2d,  coupled  imp,  lam) 
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Contours  of  Static  Temperature  (k) 
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FLUENT  6.0  (2d,  segregated,  lam) 


Contours  of  Static  Temperature  (k) 


Mar  16, 2005 
FLUENT  6.0  (2d,  segregated,  lam) 
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Contours  of  Static  Temperature  (k) 
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Contours  of  Static  Temperature  (k) 


Mar  16, 2005 
FLUENT  6.0  (2d,  segregated,  lam) 


Position  of  high-temperature  spots 


The  position  of  the  hots  spots  is  determined  by: 


•  The  relative  ratios  of  the  heat  transfer  coefficients  of 
different  tubes 

•  The  distance  of  the  smallest  tubes  from  the  slab  surface 


Next  steps 


3d  case 

Entry  lenght  analysis 

Different  heat  transfer  values  for  different 
holes  diameters 


Constructal  Workshop,  17-18  March  2005,  Rome 


Performance  Improvement  of  Tubular  Solid 
Oxide  Fuel  Cell  Stack  through  Variation  of 

its  Configuration  Design 


Michele  Call,  Chiara  Ciano,  Vittorio  Verda 
Politecnico  di  Torino  —  Dipartimento  di  Energetica 


SOFC 


Some  advantages  of  Solid  Oxide  Fuel  Cell 
high  efficiency, 
low  NOx  production, 

possible  direct  use  of  methane  by  operating  an  internal 
reforming, 

possibility  of  using  SOFC  in  combination  with  bottoming  cycle 

This  research  aims  at  studying  the  configuration  of  the  system 
integrating  cells  and  internal  reformers,  by  changing  the  cells 
geometry  and  disposition. 


From  a  single  cell  to  a  bundle 


A  model  accounting  for 
thermal  and  electrochemical 
phenomena  in  a  single  cell 

S to  calculate  the  performances 
of  the  cell  for  varyed  boundary 
conditions 

v  to  evaluate  the  effects 
produced  by  varying  geometrical 
data,  particularly  the  cell  length. 


A  network  approach  to 
compute  a  bundle 

S accounting  for  the  electrical 
connections  among  the  cells 
and  the  energy  transfers  inside 
the  bundle  and  with  the  extern 

^ to  show  how  the  system 
evolve  varying  the  shape  of  the 
bundle  and  the  number  and 
length  of  the  cells. 


System  configuration 


Three-in-parallel  by  eigth-in-seties  cell  bundles 


Electrochemical  model 


Each  cell  is  represented  as  divided  into  resistances 
connected  in  series  and  in  parallel. 


The  fuel  feeding  the  cell 
flows  for  each  layer  of  nodes 
in  branches  1,  3,  6  and  8. 

Branches  2  and  7  are  the 
electrical  connections. 

By  branches  4  and  5  the  cell 
results  electrically  isolated 
from  those  beside  it. 


Cell  interconnect 


Electrochemical  model 


Thermal  model 


Thermal  flows  are  generated  inside  the  cell  both  by  voltage 
losses  and  the  electrochemical  reaction  that  takes  place  at  the 
anode-electrolyte  interface. 


This  reaction,  consisting  in  the  oxidation  of  hydrogen  and  carbon 
monoxide,  is  the  cause  of  the  electromotive  force  E: 


Thermal  model 


Heat  is  partially  supplied  to 
the  reformers,  where  it  is 
required  for  methane  to 
produce  hydrogen  and 
carbon  monoxide. 


CHa  +  H20  ->  3 H2  +  CO 
CO  +  H2O  — ^  H2  CO 2 


In  actual  design,  reformers 
are  located  between  the 
bundles. 


Thermal  model 


A  part  of  the  thermal  flow 
released  in  the  elctrochemical 
reaction  is  recuperated  by 
heating  the  air  flux  flowing  in 
the  annular  section. 


To  reformer 


To  inlet  air 


Resistivity 

Electronic  and  ionic  resistivities  depend  on  the  local  temperature. 


Cell  simulation 


Single  cell  simulation 
(Operating  temperature  ~  900  °C) 


Current  [A] 


Boundary  conditions: 

-  DV  bundle  =  Design  DV  cell  (0.72  V)  •  Number  of  cells  in  series 

-  Fuel  temperature  (550  °C) 

-  Air  temperature  (830  °C) 

-  Heat  flow  for  internal  reforming  (30%  of  the  requirement) 

Present  design  parameters: 

-  Bundle  shape 

-  Cell  length 

Constraints: 

-  Total  active  area 

-  Anode,  cathode  and  electrolyte  thickness 


From  constructal  theory 


Actual  design 


Optimal  shape 


Resistance/design  resistance 


Electronic/ionic  resistance 


Overall  resistance  decreases  as  the  number  of  cells  in  series  increases. 


Present  and  future  work 


Geometrical  parameters  investigated/under  investigation: 

-  Cell  length  V 

-  Cell  length  and  diameter 

-  Cell  shape 


Present  work 


The  voltage  losses  decreases  as  the  cell  length  increases  (design  voltage) 


1,02 


Numerical  simulations  of  the 
fluid  &  heat  flow  in  disc¬ 
shaped  compact  heat 
echangers  with  dendritic  flow 

structure 


University  of  Roma  /-  La  Sapiensa 


Compact  heat  exchangers:  rectangular  box 
actual  geometries  multi-level  refinement 


University  of  Roma  /-  La  Sapicnza 


Constructal  compact  heat  exchanger 

committal 


A.  Be j an 
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four  levels  disc  dimensions  (mm)  1/2 


Disc  dimensions  2/2 
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Disc-shaped  constructal  compact  heat 

exchangers  model 
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Model  improvement  and  CAD  simplifications 


Discontinuities 

simplifications 


Angles  smoothed  at  the 
junctions  (15°) 


University  of  Roma  /-  La  Sapicnza 


Implications  of  the  modification 

introduced 


•  Increased  pressure 
drop  for  square 
sections 

^rect=KlX 


•  Reduced  pressure 
drops  for  smoothed 
junctions 

C=C(1-FJF,) 


(Ki=0.9) 
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Boundary  and  Initial  conditions 

Fluid:  WATER  Disc  material:  ALUMINIUM 

Outflow  Inlet  water  Temperature=300°K 

Pressure  inlet=p0+2500Pa  Inlet  hot  water  temperature  370  °K 


Pressure  outlet=p0 


Velocity  inlet=1m/s 
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Grid  Description 

Non  conformal  hybrid  mesh: 

Fluid  cells:  Hexahedra  (62900+62890)  -  solid  cells:  tetrahedra 


Level 

0 


Cells 

1539373 


Faces 

3326935 


Nodes 

450871 


Partitions 

1 
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Values  of  Reynolds  numbers  at  each 

e 


Level  0 

2548 

Level  1 

2013 

Level  2 

1605 

Level  3 

1274 

11 
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Numerical  model  details 


Solver:  segregated 
Time:  steady 
Viscous  model:  laminar 
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Velocity  profiles  -  level  0  -  ln->out 


Contours  of  Velocity  Magnitude  (m/s) 


Mar  17,  2005 
FLUENT  6.0  (3d,  segregated,  lam) 


recirculation  bubble 


‘KM  12  raw 


f-lIJifirhUi.  ■-!.  li^i 


7.89e-01 

7.17e-01 

6.46e-01 

5.74e-01 

5.02e-01 

4.30e-01 

3.59e-01 

2.87e-01 

2.15e-01 

1.43e-01 

7.17e-02 

O.OOe+OO 


University  of  Roma  /-  La  Sapienza 


j  t  fc^|Bf.wpp|ig  nj  1 1 1 1 1 1  f*ttM J. WHWfMWVW 


1fPrtf¥MfXKWWWi 


Temperature  profiles  -  level  0  -  ln->out 


amrmcml 


3JH"TZ 
: :  n-3  :■ 
]»M^ 

UJr-JZ 


University  of  Roma  /-  La  Sapienza 


Zoom  on  the  recirculation  bubble 

(temperature) 
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Temperature  profiles  -  level  0  -  ln->out 
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Temperature  profiles  -  out->in 
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Temperature  profiles  at  the  hot  fluid  inlet 
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Contours  of  Static  Temperature  (k) 
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Contours  of  Static  Temperature  (k) 

Mar  12,  2005 
FLUENT  6.0  (3d,  segregated,  lam) 

performance 
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Area-Weighted  Average 


Static  Temperature  (k) 


outflowl  321.80 

outflow2  321.67 

outflow3  320.00 

outflow4  320.32 


Area-Weighted  Average 
Static  Temperature  (k) 


central  outflow  350.30 


Area-Weighted  Average 

Static  Pressure  (pascal) 


outflowl  -836.25 

outflow2  -1018.8 

outflow3  -978.87 

outflow4  -748.10 
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Proposal  for  an  improved  model: 

Tubes  slightly  curved  to  obtain  symmetrical 
junctions 

Exploit  to  the  upmost  the  boundary  layer 
development  along  each  tube 

Splits  optimized  on  the  basis  of  pressure  drops 


University  of  R  oma  b  la  Sapienza 
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Solution  Hystograms 


Dendritic  pattern 
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Performance  Graphs 
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Pressure  drop  on  a  diverging  section 


Simulation  and  measurements  of  turbulentheat  transfer  in  a  channel  with  a  surface 
lounted  heated  block  (chen-  wang),  heat  and  mass  transfer,  31 ,  1996,  463-473- 
Springer  Verlag 
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Energy  usage  in  Aircraft  Equipment  Systems 

A  potential  future,  from  the  EC  Power  Optimised  Aircraft  project 


Workshop  on  Constructual  Theory 
Rome,  17-18  March  2005 


Lester  Faleiro,  Jacques  Herzog 
Liebherr-Aerospace,  Germany 


Project  Goals 


The  target  of  POA  is  to  validate,  at  aircraft  level  and  both  qualitatively  and 
quantitatively,  the  ability  of  next  generation  aircraft  equipment  systems  to  enable 
the  reduction  in  consumption  of  non-propulsive  power 

►  Drivers 

□  Safety  Standards 

►  Objectives 

□  Reduction  of  peak  non-propulsive  power  by  25% 

□  Reduction  of  total  non-propulsive  power 

□  Reduction  of  fuel  consumption  by  5% 

□  Reduction  of  total  equipment  weight 

►  Constraints 

□  Maintenance  Costs 

□  Equipment  production  costs 

□  Reliability 


LIEBHERR-AEROSPACE  2005 


DO-05-TFT-0302  /  17.03.2005 
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Conventional  Aircraft  Architecture 


Energy  usage  in  Aircraft  Equipment  Systems 


LIEBHERR-AEROSPACE  2005 
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Optimised  Aircraft  Architecture? 


r 


Primary 
Controls 


Secondary 
Controls 


Environmental  ^ 
Control 


Commercial 
Loads 


_*i*' 

i**  ^ 


Cabin  Expansion 
generator 


Landing 
Gear 


No  Gearbox 


Local  Compressor 


Reduced  Engine  Bleed 


Local  Hydraulic  source 


More  Electrical  Power 


Energy  usage  in  Aircraft  Equipment  Systems 


LIEBHERR-AEROSPACE  2005 
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Electric  Actuators 


Fan  Shaft  Generator 


•  ~  150kW  main  Power  Generation  at  Idle  and  Above 

•Emergency  Power  Generation  from  Windmilling  Fan 
to  Airframe  (engine-out)  or  HP  Machine  for  Assisted 
Windmill  Relight. 


•High  Pressure 
Starter/Generator 


•  -  200kW  Motor  for  Engine  Starting. 

•Will  Generate  Power  after  Engine  Start 

•LP  to  HP  Power  Transfer  may  Improve  Fuel 
Burn  and  Reduce  Thrust  During  Descent. 

THALES 


Goodrich. 


DC  Power  Bus  on  Engine 

•Simplified  Airframe/Engine  Interface 

•Each  Machine  will  have  a  Power  Electronic 
Drive 

•Each  Drive  will  Appear  as  a  Node  in  a 
Distributed  Control  System 


Active  Magnetic  Bearing 

•Investigate  Potential  for  Removal  of  Oil 
System 

•Monitoring  of  Shaft  Rotordynamics 


Electric  Oil  Pump/Scavenge  Syste 

•Optimise  Oil  Flow  to  Bearings  over  Engine  Cycle 

Techspace  Aero 

snecma  group 


Electric  Fuel  Metering  Unit 

100kW  Motor 
•Simplified  Fuel  System 
•Lower  Heat  Input  to  Fuel 


Electric  Oil  Breather  Model 

zz^*)  Hispano-Suiza  QA/$ 

groupe  snecma  Y,  Z  Y 


AVIONICS  ELECTRICAL  SYSTEM 


groupe  snecma 


Actuation  Systems 

The  Objective  is  to  achieve  lower  life  cycle  costs,  through 
power  optimisation,  reduced  weight  and  maintenance  costs 


Goodrich 


Electro-Hydraulic  Actuation  (EHA) 

Derivation  of  Standards 
Hardware  and  model  verification 


2.5  kW  EHA 


Nacelle  Systems 


More  Electrical  Thrust  Reversal 

Alternative  to  pneumatic  and  hydraulic  thrust 
reversal  systems 

Mechanical  and  data  bus  synchronisation 


^pp)  Hurel-Hispano 


s necma  group 


Goodrich 


Snecraa  IVfoteurs 


grmpc  SHcontj 


Electro-Back-up  Hydraulic 
Actuation  (EBHA) 

Novel  and  versatile  hybrid  actuation 


Trimmable  Horizontal 
Stabiliser  Actuation  (THSA) 

Proof  of  concept 

More  electrical  actuation  with  innovative 
mechanical  technologies 

Goodrich  insu 


§ 

Wide-body  aileron  actuator  ~2  kW 


Typical  wide-body  stabiliser  actuator 


A 


Spoiler  actuator  ~25  kW 


Electro-Mechanical 
Actuation  (EMA) 

Distributed  High-Lift  systems 


Comparison  of  hinge  line 
versus  rotary  technologies 


Goodrich 


Landing  Gear  Systems 

Landing  Gear  system  integration 
Decentralised  actuation  for  Nose  Wheel 
Airbus  A300/A310  More  Electrical  Actuation  for  Main  Gear 
More  Electrical  Wheel  Braking 


Messier-Dowty 

SHt’fjjfdZ  group 


Labinal 

gr&upe  snecma 


Pneumatic  Systems 


The  Objective  is  to  reduce  and  optimise  the  effect  of  bleed  air  off-take  on  power 
usage,  as  this  is  a  large  consumer  of  non-propulsive  power.  Bleed  Air  Off-Takes 
are  mainly  used  for  the  Environment  Control  System  and  for  Wing  Ice 

Protection. 


Conditioning  System  (ACS) 

To  increase  the  efficiency  of  the  ACS 

The  main  innovation  is  the_combination  of  a  Vapour  Cycle  (containing 
an  environmentally  neutral  fluid)  with  an  electrical  driven  Air  Cycle.  A 
variable  speed  motor  for  the  re-circulation  fan  and  a  Cabin  Energy 
Recovery  Device  will  be  used. 

The  main  outputs  are  a_Model  of  a  complete  ACS  and  the  Test  of  a 
Hybrid  ACS  (Vapour  +Air  Cycle) 

TORiNCMTALV 


Re-circulation  Fan 


Motorised  Air  Cycle  machine 


Wing  Ice  Protection  (WIP) 


To  Increase  the  efficiency  of  WIP  Systems 

The  main  innovation  is  the  use  of  ultrasonic  surface  ice  sensors  (Anti  Ice  on 
demand)  and  hybrid  wing  heating  (electrical  and  hot  air) 

The  main  outputs  are  Models  of  WIP  Systems  and  Test  of  an  innovative  WIP 
System 

Wing  heat  distribution  using:  Ultrasonic  sensors,  Electro-thermal  devices, 
On  demand  active  intelligence  control  and  Monitoring  of  unprotected  surfaces 


Ice  detection  sensor 


Fuel  Cells  (FC) 


To  Increase  the  efficiency  of  electrical  power  generation  for 
pneumatic  systems 

The  main  innovation  is  the  validation  of  a  Solid  Oxide  Fuel  Cell 
(SOFC)  with  its  reformer  for  use  with  kerosene 

The  main  outputs  are  a  Model  of  FC  System  and  the  Test  of  a  5  kW 
Fuel  Cell  System 


Fraunhofer  |n5tjtut 

So  I  a  re  Energiesysteme 


B  Alenia 

AERONAUTICS 


%,\4s»SK!R> 


SOFC  Principle 


Project  Process  so  far 


►  The  A330-300  was  chosen  as  a  reference  for  all  comparison 

►  More  Electrical  systems  chosen  as  technology  solutions  for  POA 

□  Perception  of  more  efficiency  at  the  system  level 

□  Higher  synergy  when  integrated  into  a  more  electrical  aircraft  (MEA) 

►  Systems  combined  into  various  POA  (mostly  MEA)  architectures 

□  Each  system  has  been  designed  and  many  are  being  manufactured 

►  Note:  This  is  NOT  optimisation  -  this  is  a  pre-selection  of  a  few  possible 
combinations  that  the  industrial  partners  together  feel  have  the  greatest 
potential  for  application  in  the  near  to  medium  future 
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Energy  related  trends  at  the  Systems  Level 


►  MESys  tend  to  waste  less  energy 

□  No  hydraulic  leakage 

□  No  large  losses  from  engine  bleed 

□  Power  on  demand 

►  Many  technologies  can  be  demonstrated  today 

□  MEE  technologies 

□  High  density  power  generation 

□  Reliable  electrical  distribution  architectures 

□  Better  aeronautical  cabling 

□  Electro-hydrostatic  and  electro-mechanical  flight  control  actuation 

□  Electrical  braking 

□  Electrical  Environmental  Control  System 
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One  challenge  for  Constructal  Theory  -  The  ECS 


Energy  usage  in  Aircraft  Equipment  Systems 
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Tasks  of  the  Environmental  Control  System  (ECS) 


►  Pressurisation  (at  altitude  Pressure  not  lower  than  752  hPa) 

►  Ventilation  (Fresh  flow  per  passenger  0,55  Ib/min) 

►  Temperature  control  (Cabin  temperature  around  24°C) 

□  Engine  Bleed  Air  and  Pre-cooling 
with  Engine  fan  air 

□  Bleed  Air  with  Ram  Air 

□  Turbomachines  and  Water 
separation  loops 

□  Mixing  chamber  and  distribution 

Fig.l 
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Principle  of  the  current  ECS  on  large  civil  Aircraft 


► 

► 


► 


2  bleed  ports:  LP  port  at  high  engine 
power  and  HP  port  at  low  engine  port. 

Pre-cooler  limits  the  temperature  to 
200°C  using  Fan  air 

PRV  regulates  the  pressure  delivery 


►  Flow  controlled  by  FCV 

►  Temperature  controlled  by  Ram  air  flow 
(RAI)  and  Bypass  flow  (BPV) 

►  Water  extraction  realised  by  high 
pressure  water  separation  loop  (HPWS) 
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Energy  Impacts  of  the  current  ECS  on  a  large  civil  Aircraft 


►  Impact  on  the  Engine  -  High  Specific  Fuel  Consumption 

►  Losses  -  high  temperature  and  pressure  air  has  to  be  cooled  and  de-pressurised 
before  it  is  useable 

□  Weight  of  components  included  in  order  to  induce  energy  losses 

►  Drag  for  the  cooling  of  the  Ram  Heat  Exchangers 

►  And  many  more.... 
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Power  losses  due  to  the  ECS  on  current  Aircraft 


►  Currently  the  Engines  provide  (pneumatic)  power  to  the  ECS.  This  supplied  power  is  strongly  dependent  on  the  flight  phase. 

►  Two  bleed  ports  needed  to  supply  the  required  power  and  to  limit  the  impact  of  Bleed  Air  on  Engine  fuel  consumption 

►  Nevertheless  a  lot  of  losses  are  induced  by  the  divergence  of  needs  and  supply 
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ONE  solution  to  energy  “optimisation”:  the  electrical  ECS  (Bleed-less) 


mechanical  energy  electrical  energy 


Generators 


ECS 


pneumatic  energy 


Cabin 


►  Engines  produce  mechanical  (kinetic)  energy  from  fuel 

►  Generators  convert  the  mechanical  energy  into  electrical  energy 

►  ECS  use  the  electrical  energy  to  perform  its  tasks: 

□  Pressurisation 

□  Ventilation 

□  Temperature  control 
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Adaptability  of  the  electric  ECS 


*  The  power  consumption  of  the  electrical  ECS  is  independent  of  the  engine  speed  (no 
dissipation  of  excessive  engine  power) 

►  Adaptation  of  the  ECS  to  the  cabin  needs 

□  Fresh  flow  can  be  adapted  to  the  number  of  Cabin  occupants 


□  Cooling  power  can  be  adapted  to  the  heat  loads 
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Impact  on  engine  of  pneumatic  and  electrical  power  source 


Mechanical  Off-takes  (to  produce  electrical  power)  and  Pneumatic  Off-takes  have  different 
impacts  on  the  Specific  Fuel  Consumption  of  an  engine 

□  to  cover  the  ECS  needs,  the  impact  of  electrical  power  on  the  engine  specific  fuel  consumption  is  not 
so  strong  as  the  impact  of  the  pneumatic  off-takes 


Potential  to  reduce  the  Fuel 
consumption  (energy)  by 
using  electrical  power -this 
pseudo-optimisation  is  what 
goes  on  at  the  equipment 
systems  level 


Thrust 


Fig.  11 
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A  role  for  Constructal  Theory? 


►  A  Bejan  has  published  2  papers  to  address  the  issues  in  an  aircraft  ECS 

□  The  method  is  structured  and  makes  sense,  but  the  problem  to  solve  is 
complex  and  multi-faceted. 

•  It  cannot  be  solved  by  only  thermodynamic  considerations,  as  the  systems 
involved  are  multi-physics. 

•  The  inputs  required  to  address  the  problem  fully  are  difficult  to  obtain,  as 
many  organisational  entities  may  be  involved  (aircraft,  engine  and  systems 
manufacturers) 

•  It  is  difficult  to  solve  problems  for  one  component  without  addressing  the 
impacts  on  the  other  components 

•  It  is  difficult  to  solve  problems  for  the  ECS  without  addressing  the  impacts 
on  engine  and  rest  of  the  aircraft  (drag,  etc.) 

►  Could  Constructal  Theory  produce  a  technological  solution  such  as 
bleedless  ECS,  or  is  something  else  needed  to  complement  it,  such  as  a 
palette  of  potential  technological  solutions  and  a  palette  of  constraints? 
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Further  Challenges  for  Energy  Based  Methods 


►  Constraints  are  not  easy  to  include  in  such  methods 

□  Weight  limits 

□  Volume  limits 

□  Design  limitations  on  the  performance  and  function  of  Aircraft 
Equipment  Systems 

□  Safety  considerations 

►  Every  divergence  from  the  ideal  model  can  reduce  the  degree  of 
freedom  of  the  method 

□  The  idea  of  creating  an  ideal  energy  based  model  and  then 
systematically  altering  the  model  to  include  design  constraints  is  good, 
but  the  complexity  of  such  a  model  is  enormous 
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Energy  Based  Methods  are  probably  needed 


►  For  a  really  Power  Optimised  Aircraft  (not  just  MEA),  thinking  must  start  to 
occur  at  the  level  of  functions,  rather  than  systems 

□  Design  Process  Re-engineering  for  equipment  systems  integration  may  be 
necessary  -  energy  based  methods  might  be  useful 

□  e.g.  ECS  is  currently  used  to  provide  pressurisation  and  cooling.  In  POA, 
various  architectures  that  provided  these  two  functions  in  various  numbers  of 
systems  and  sub-systems  were  examined 

►  Decreasing  Systems  Autonomy  -  energy  is  the  common  denominator 

□  More  Electrical  Engine  (MEE)  and  MEA  go  hand-in-hand 

□  Availability  of  MEA  engine  thrust  is  connected  to  generation  system  reliability 
through  the  use  of  electrical  fuel  pumps 

□  This  trend  on  interdependency  is  most  visible  on  Boeing  B787:  bleedless 
engine/electrical  ECS/electrical  generation  and  distribution  -  can  energy  based 
methods,  such  as  Constructual  Theory,  come  up  with  such  solutions? 
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New  Challenges 


►  Cooling  in  a  MEA 

►  Thermal  management  of  Power  Electronics  is  a  challenge 

►  MOET  Project  -  European  University  to  look  at  this? 
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Thank  you! 
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On  the  synthesis  of  complex  energy  systems  with  several  heat  transfer  interactions 
by  Andrea  Lazzaretto 


Synthesis,  Analysis  and  Optimization 


(Bejan,  Moran,  Tsatsaronis,  1996  'Thermal  Design  and  Optimization",  Wiley  ed.) 


•  SYNTHESIS  is  concerned  with  putting  separate  elements  into  a  whole,  i.e. 
the  particular  equipment  items  making  up  the  overall  thermal  system  and 
their  interconnections  are  specified 

•  The  objectives  of  the  ANALYSIS  and  OPTIMIZATION  are  to  identify  the 
preferred  configuration  among  those  synthesized. 
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Synthesis,  Analysis  and  Optimization 

•  ANALYSIS  generally  entails  thermal  analysis  (mass,  energy  and  exergy 
analysis),  costing  and  sizing  of  equipment  on  at  least  a  preliminary  basis 

•  OPTI  Ml  ZATI  ON  can  take  two  general  forms: 

1  )  PARAMETER  and  2)  STRUCTURAL  optimization 

1)  pressures,  temperatures,  mass  flow  rates  and  chemical  composition  at 
various  points  and/or  other  system  key  variables  are  determined,  at  least 

approximately,  with  the  aim  of  satisfying  some  desired  objectives 

2)  the  equipment  inventory  and/or  interconnections  among  components  are 
altered  iteratively  to  achieve  a  superior  design 
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Synthesis,  Analysis  and  Optimization 

Concept  Generation  and  ^ 

Preliminary  Screening _ — J 


Detailed  Design 
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Other  criteria  to  improve  thermal  system  structure 

•  Definition  of  a  superstructure  including  all  feasible  structures  and 
optimization  of  design  parameters  to  "cut"  useless  branches  of  the 
superstructure 

•  Criterion  applied  here  is  based  on  a  "basic  plant  configuration" 

which  is  defined  by  subdividing  the  plant  into  two  zones: 

-  One  including  the  " basic  components  that  are  always  included  and 
connected  in  the  same  way  in  the  plant  structure; 

-  The  heat  transfer  zone,  considered  as  a  " black  boS  in  which  heat  is 

transferred  between  hot  and  cold  thermal  flows  (made  available  from  the 
components  in  the  first  zone)  independently  of  the  number,  type  and 
interconnections  among  heat  exchangers,  which  are  not  defined. 
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Optimization  of  the  basic  plant  configuration 

•  General  problem:  optimal  synthesis  of  the  components,  i.e.  how  to 
interconnect  them  in  order  to  reach  the  optimal  value  of  the  objective 
function  (here  the  exergetic  efficiency) 

•  Idea:  fix  a  part  of  the  structure  (basic  components  and  their 
interconnections)  and  keep  the  definition  of  the  heat  exchanger  structure 
independent  of  the  structure  of  the  overall  system 

•  Suggested  solution:  parameter  optimization  of  the  basic  plant 
configuration  +  construction  of  the  HEN  which  fulfill  the  optimal 
thermodynamic  conditions  at  the  black-box  boundaries  and  guarantees 
the  feasibility  of  the  heat  transfer 

•  Final  solution:  find  a  criterion  to  generate  all  possible  basic  plant 
configurations  and  then  choose  and  optimize  the  best  one  according 
to  the  objective  of  the  analysis 
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HAT  cycle  "  Basic  Plant  Configuration  " 

Lazzaretto,  Segato,  2001,  Thermodynamic  Optimization  of  the  HAT  cycle  Plant,  J  for  Eng.  of  Gas  Turbine  and  Power 
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HAT  cycle  "  Basic  Plant  Configuration  " 
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Heat  transfer  interactions 
within  the  basic  plant  configuration 


On  the  synthesis  of  complex  energy  systems  with  several  heat  transfer  interactions 
by  Andrea  Lazzaretto 


Heat  transfer  section  (black-box) 


A  A  A  |  A  A  A  A 


Hot  side 


Cold  side 
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HAT  cycle  "  Basic  Plant  Configuration  " 

•  Basic  components:  compressors  (cl,  c2),  combustion  chamber  (cc), 
turbine  (t),  saturator  (s) 

•  Order: 

air/gas  loop:  (cl,  c2),  (cc),  (t) 
water  loop  (s),  (cc),  (t) 

•  Temperature  links  between  basic  components  are  “virtually"  cut, 
i.e.  the  temperature  at  the  inlet  of  a  component  is  left  independent 
of  the  temperature  at  the  outlet  of  the  preceding  component 

(with  the  exceptions  of  the  links  between  hot  source  (cc)  and 
turbine  (t)  and  between  cold  source  (ambient)  and  compressor  (cl) 
where  additional  heat  transfer  is  not  meaningful) 
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Structure  of  complex  energy  systems 


•  One  or  more  thermodynamic  cycles  are  involved 

•  Each  cycle  involves  a  loop  of  one  or  more  materials 

•  I  n  each  cycle  the  order  of  the  basic  components  is  defined 

•  Additional  basic  components  can  be  added  to  improve  the  system 

performance  (e.g.  compressor  and/or  turbine  stages,  combustion 
stages) 
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Parameter  optimization  (thermodynamic) 
of  the  basic  plant  configuration 

•  Objective  function:  total  plant  exergetic  efficiency 

•  Decision  variables  are  the  degrees  of  freedom  needed  to  set  the 

thermodynamic  conditions  at  the  inlet/outlet  of  the  basic  plant 
components  -> 

once  their  values  are  fixed,  inlet  and  outlet  temperatures  and  mass 
flow  rates  at  the  black-box  boundaries  are  fixed 

•  Optimization  problem  is  subject  to  the  heat  transfer  feasibility 
constraint,  simply  expressed  by  the  condition  that  the  minimum 
temperature  difference  between  composite  curves  should  be 
greater  than  a  minimum  AT 
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Composite  Curves 

•  Thermal  capacities  of  flows  which  transfer  heat  in  the  same 
direction  (hot  to  cold  or  vice  versa)  and  in  the  same  temperature 
interval  are  summed  up. 

•  A  Hot  and  Cold  Composite  Curves  (HCC,  CCC)  are  so  drawn  in  a 
heat-temperature  diagram,  showing  the  sequence  of  the  heat 
transferred  at  increasing  temperature  values  according  to  the 
thermal  capacities  calculated  for  hot  and  cold  flows,  respectively 

•  No  external  hot  sources  exist  ->  HCC  and  CCC  have  the  same 
ending  abscissa 

•  The  black-box  may  release  heat  to  the  ambient  ->  the  projection  of 
the  HCC  on  the  abscissa  is  longer  than  that  of  the  CCC 


Temperature  [°C] 
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Composite  Curves 
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Heat  transfer  in  the  black-box 


•  In  the  search  for  the  optimum  total  plant  exergetic  efficiency,  the 
optimization  algorithm  finds  implicitly  the  condition  of  most 
efficient  heat  transfer  in  the  black-box  (optimum  matching  between 
hot  and  cold  thermal  flows) 

->  optimum  matching  is  meant  as  the  set  of  inlet  and  outlet 
temperatures  and  mass  flow  rates  of  the  hot  and  cold  streams 
deriving  from  the  independent  variables  set  which  minimizes  the 
distance  between  the  CC 
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Temperature  Profiles 
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HAT  cycle  heat  exchanger  network  ( 1) 


PINCH  POINT 
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HAT  cycle  heat  exchanger  network  (2) 
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Heat  exchanger  network  (3) 


PINCH  POINT 
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Heat  exchanger  network  (4) 


PINCH  POINT 
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HAT  cycle  plant 
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Building  the  heat  exchanger  network  (HEN) 


Need  of  fulfilling  black-box  optimal  boundary  conditions  (corresponding  to 
achieving  the  maximum  total  plant  exergetic  efficiency)  and  minimum  pinch 
point  (ATpp) 

HEN  can  be  built  in  two  different  ways: 

1)  No  mixing  at  different  temperatures 
->  same  composite  curves; 

2)  Mixing  at  different  temperatures 

->  different  composite  curves  (ATPP  may  or  not  be  fulfilled) 


Other  solutions  in  which  ATPP  is  respected  but  boundary  conditions  are  not 
fulfilled  ->  lower  exergetic  efficiency  -  different  composite  curves 
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Basic  Plant  Configuration 
Regenerative  gas  turbine  cycle 


compressor 
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Air  preheating  -  composite  curves 
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CRCC  "  Basic  Plant  Configuration  " 
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CRCC  heat  transfer  section 


Fig.2  Black-box  representing  the  heat  transfer  section  in  the 

basic  plant  configuration 
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I  nteractions  between  HRSG  and  black-box 
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Design  improvement  procedure  for  the 
"CRCC  Basic  Plant  Configuration" 

•  Problem:  high  number  of  free  variables  in  the  design  improvement 
procedure: 

1)  decision  variables  in  the  components  outside  the  heat 
transfer  section 

2)  variables  at  the  boundaries  of  the  heat  transfer  section 
(inlet/outlet  temperatures  and  mass  flow  rates  of  the  thermal 
flows  in  the  black-box) 

->  need  of  reduction  of  the  degrees  of  freedom  at  the  heat 
transfer  section  boundaries 

->  Because  of  the  system  complexity,  system  is  simulated  into 
parts  (sections),  no  optimization  algorithms  are  used 
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Design  improvement  procedure  for  the 
"CRCC  Basic  Plant  Configuration" 


•  Preliminary  thermodynamic  analyses  are  performed  to  understand 
whether  some  of  the  decision  variables  can  be  fixed  as  parameters  or 
show  predictable  trends 

•  Series  of  simulation  runs  of  the  total  system  are  carried  out  by  varying 
the  remaining  decision  variables  one  at  a  time 

•  Basic  importance  of  a  correct  choice  of  the  reformer  operating  variables 


Efficiency  Unconverted  Methane  Ratio  [%] 
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Preliminary  thermodynamic  analysis 


Fraction  of  exhaust  gases  sent  to  the  reformer  CC 


HTS  Temperature  [°C] 


LTS  Temperature  [°C] 
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I  ndependent  variables  of  the 
“basic  plant  configuration" 


VARIABILE 

TARGET 

Reformer  Pressure 

PREF 

Optimize  (10  -  22-  30  bar) 

Gas  Temperature  at  the  reformer 
outlet 

Steam  /  Natural  Gas  Ratio 

Tref 

SGR 

Tref  and  SGR  are  functions  of  pRtF 
and  they  are  related  if  methane 
conversion  ratio  is  fixed  (95%) 

HT  S  T  emp  erature 

Thts 

Optimize  (range:  350-475°C) 

LT  S  T  emp  erature 

Tlts 

Optimize  (range:  180-250°C) 

Water  temperature  at  the  saturator 
inlet 

T  in, SAT 

Optimize 

Water  mass  flow  rate  at  the 

saturator  inlet 

min,SAT 

Fix  at  constant  value  (mole  frac.  H2) 

Exhaust  gas  mass  flow  rate  to 
reformer 

mPCC 

Optimize  (range:  35-100%  GT  outlei 

Reformer  reactive  mixture 

temperature 

Tr 

Optimize  (450  -  500  -  565°C) 

Exhaust  gas  temperature  at  the 
reformer  outlet 

TfUMI,REF 

Optimize  (range:  480  -  540°C) 

Preheating  of  SYNGAS 

TpRE, SYNGAS 

Minimize 

Preheating  of  natural  gas 

TpRE,NG 

Optimize  (range:  500-800°C) 

IP  steam  mass  flow  rate 
superheated  with  heat  recovered  in 
the  black-box 

mip 

Optimize  (range:  10-20  kg's) 

HP  water  mass  flow  rate  heated, 
evaporated  and  superheated  with 
heat  recovered  in  the  black-box 

mHp 

Maximize 
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Heat  transfer  maximization 

•  I  n  order  to  have  the  maximum  feasible  heat  transfer  in  the  “black¬ 
box",  a  Fortran  routine  is  included  in  the  plant  simulation  model  to 
set  the  last  of  the  decision  variables  (mHP)  while  keeping  the  other 
variables  as  constant  parameters  varied  in  selected  ranges 

As  in  previous  examples: 

•  the  heat  transfer  feasibility  constraint  applies,  i.e.  the  minimum 
temperature  difference  between  composite  curves  should  be 
greater  than  a  minimum  AT  (it  is  directly  applied  to  solve  for  mHP) 

•  no  external  hot  sources  exist  ->  HCC  and  CCC  have  the  same 
ending  abscissa 

•  the  black-box  may  release  heat  to  the  ambient  ->  the  projection  of 
the  HCC  on  the  abscissa  is  longer  than  that  of  the  CCC 
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Heat  transfer  maximization 

•  The  only  HP  water  mass  flow  rate  re- circulated  outside  the  HRSG  (mHP) 
is  considered  as  decision  variable  in  the  black-box  sub-optimization  (heat 
transfer  maximization)  while  keeping  the  other  variables  as  constant 
parameters  (varied  in  selected  ranges) 


•  Each  mHP  value  corresponds  to  a  specific  shape  and  position  of  the  Cold 
Composite  Curve.  1000 
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Fig.  9  Cold  composite  curves  obtained  using  values  of  mHP  lower  than  (CCC1),  greater  than 
(CCC2)  and  equal  to  (CCC)  the  maximum  one. 
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Decisions  about  thermal  flows 
to  be  included  in  the  black-box 

•  In  general,  to  have  maximum  freedom  in  the  definition  of  the 
system  structure  and  parameters,  all  thermal  flows  should  be 
included  in  the  heat  transfer  section  (black-box) 

Thermal  flows  not  “directly"  included  in  the  black-box  of  the  CRCC 
plant: 

•  structure  of  the  HRSG  is  considered  to  be  fixed  (is  part  of  the  basic 
components  of  the  system),  but  HP  and  I P  steam  can  be  produced 
also  outside  the  HRSG  using  heat  available  from  the  black-box  (p 
and  T  levels  decided  in  advance) 

•  p  and  T  levels  of  thermal  flows  taken  from  the  black-box  and  used 
in  the  Section  for  Absorption/Stripping/Compression  of  the  C02 
decided  in  advance 
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I  nfl uence  of  variations  of  the  design 
variables  on  the  composite  curves 


0  100000  200000  300000  400000  500000  600000  700000  80000 

AH  [kW] 

Fig.  12  "Composite  Curves"  at  maximum  total  plant  efficiency  for  three  different  thermal 
profiles  (numbered  according  to  the  order  in  Figure  10)  of  the  reformer  reactor  operating 
at  22  bar. 


[3-11 
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I  nfluence  of  variations  of  the  design 
variables  on  the  black-box  composite  curves 


AH  [kW] 


Fig.  11  "Composite  Curves"  for  two  values  of  the  water  temperature  at  the  saturator  inlet 
and  reformer  pressure  fixed  at  22  bar. 
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Priority  order  in  the 
variables  in  the  First 


variation  of  the  decision 
Series  of  simulation  runs 


Outlet  temperature,  SGR  and 
thermal  profile  for  the  3  reformer 
sections  (95  %  CH4  conversion): 

(a)  840°C/ 3.5/ 840-720-570 

(b)  800°C/ 5/ 800-680-530 

(c)  775°C/ 6/ 775-655-505 

HTS  temperature:  400°C 
LTS  temperature:  181°C 


Water  temperature  at  the 
inlet  of  saturator:  it  is 
function  of  the  saturation 
temperature  at  the  defined 
pressure  ->  tmin=140°C 

140, 160, 180, 200, 220°C 


Exhaust  gases  mass  flow  rate 
to  reformer  combustor 
35-75% 

Fixed  at  50% 

Reactive  mixture  temperature 
450°C,  500°C,  565°C 
Fixed  at  450°C 

Exhaust  gases  temperature  at 
the  reformer  outlet 
480°C,  500°C,  520°Q  540°C 
Fixed  at  480°C 

SYNREF  temperature 
500°C,  600°C,  700°Q  800°C 
Fixed  at  500°C 

Natural  gas  preheating  temp. 
200°C,  300°C,  400°C,  500°C 


Outlet  temperature,  SGR  and 
thermal  profile  for  the  3  reformer 
sections  (95%  CH4  conversion): 

(a)  920°C/ 3.5  /  920-800-650 

(b)  900°C/  4/  900-780-630 

(c)  850°C/ 6  /  850-730-580 

HTS  temperature:  400PC 
LTS  temperature:  181°C 


Water  temperature  at  the 
inlet  of  saturator:  it  is 
function  of  the  saturation 
temperature  at  the  defined 
pressure  ->  tmin=175°C. 

175, 195, 215, 235°C 


Exhaust  gases  mass  flow  rate 
to  reformer  combustor 
35-75% 

Fixed  at  50% 

Reactive  mixture  temperature 
450°Q  500°Q  565°C 
Fixed  at  450°C 

Exhaust  gases  temperature  at 
the  reformer  outlet 
480°C,  500°C,  520°C,  540°C 
Fixed  at  480°C 

SYNREF  temperature 
500°C,  600°Q  700°Q  800°C 
Fixed  at  500°C 

Natural  gas  preheating  temp. 
200°Q  300°Q  400°C,  500°C 


Outlet  temperature,  SGR  and 
thermal  profile  for  the  3  reformer 
sections  (95%  CH4  conversion): 

(a)  950°C/  3.5  /  950-830-680 

(b)  900°C/ 5  /  900-780-630 


HTS  temperature:  400°C 
LTS  temperature:  181°C 


Water  temperature  at  the 
inlet  of  saturator:  it  is 
function  of  the  saturation 
temperature  at  the  defined 
pressure  -^tmin=185°C. 


185, 205, 225, 245, 265°C 


Exhaust  gases  mass  flow  rate 
to  reformer  combustor 
35-75% 

Fixed  at  50% 

Reactive  mixture  temperature 
450°C,  500°C,  565°C 
Fixed  at  450°C 

Exhaust  gases  temperature  at 
the  reformer  outlet 
480°C,  500°C,  520°C,  540°C 
Fixed  at  480°C 

SYNREF  temperature 
500°C,  600°C,  700°C,  800°C 
Fixed  at  500°C 

Natural  gas  preheating  temp. 
200°C,  300°C,  400°C,  500°C 
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Design  improvement  procedure: 
First  series  of  simulation  runs 


Design  variables  are  varied  in  the  ranges  of  Figure  10,  and 


•  the  intermediate  pressure  mass  flow  rate  re-circulated  outside  the  HRSG  (mlp 
=  20  kg/s),  and 


•  the  H2  molar  fraction  at  the  saturator  outlet  equal  to  65% 

The  high  pressure  mass  flow  rate  re-circulated  outside  the  HRSG  calculated  by 
the  Fortran  routine  in  order  to  maximize  the  heat  transfer  in  the  black  box  is 
mHP  =  115  kg/s  and  the  maximum  thermal  efficiency  of  the  total  plant  is 
50.96% 
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I  ndications  for  design  improvement 

from  the  composite  curves 
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Fig.  13  Possible  actions  on  the  cold  composite  curve  (CCC)  for  improving  the  performances 
of  the  best  plant  configuration  obtained  in  the  first  series  of  simulation  runs  (reformer 
pressure  equal  to  30  bar). 


On  the  synthesis  of  complex  energy  systems  with  several  heat  transfer  interactions 
by  Andrea  Lazzaretto 


Second  Series  of  simulation  runs 

•  I  n  the  second  series: 

-  The  intermediate  pressure  mass  flow  rate  (mlp)  re-circulated  outside  the 
HRSG  is  decreased  from  20  to  10  kg/s; 

-  The  H2  molar  fraction  at  the  saturator  outlet  is  increased  from  65%  to 
70%. 


Among  the  variables  indicated  in  Figure  10: 

-  The  preheating  temperature  of  the  syngas  to  the  reformer  before 
expansion  is  set  to  250° C,  instead  of  500° C 

-  The  temperature  of  the  exhaust  gases  from  the  reformer  is  set  at  400°C, 
instead  of  480° C 

•  Using  these  values,  the  high  pressure  mass  flow  rate  re-circulated  outside 
the  HRSG  obtained  by  the  Fortran  routine  to  maximize  the  heat 

transfer  in  the  black  box  is  mHP  =  120  kg/s  instead  of  115  kg/s  and 

the  maximum  thermal  efficiency  of  the  total  plant  is  52.32%. 
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Results  of  the 

design  improvement  procedure 


Series 

of 

Simul. 

Runs 

Pref 

[bar] 

SGR 

Temperature 
Profile  [°C] 

Tsat 

[°C] 

PXOt 

[MW] 

Pspec  [kJ/kg] 

C02  Spec. 
Emissions 
[kgC02/kWh] 

Tit  [%] 

Reference  Simple  Combined  Cycle  Configuration 

/ 

/ 

/ 

/ 

396 

717.4 

0.363 

57.76 

Optimization  Steps 

3.5 

570-720-840 

200 

392 

770.71 

0.061 

47.25 

10 

5 

530-680-800 

200 

383 

763.3 

0.057 

46.16 

6 

505-655-775 

200 

375 

751 

0.058 

45.22 

I 

3.5 

650-800-920 

215 

405 

870.17 

0.055 

48.98 

22 

4 

630-780-900 

215 

404 

832.4 

0.054 

48.82 

6 

580-730-850 

215 

395 

833.6 

0.051 

47.65 

30 

3.5 

680-830-950 

245 

410 

847.54 

0.058 

49.37 

5 

630-780-900 

245 

423 

893.7 

0.0514 

50.96 

Final  Optimized  Configuration 

II 

30 

5 

630-780-900 

245 

434.16 

900 

0.05 

52.32 

Tab.  4  Results  obtained  by  the  first  (I )  and  second  (II)  series  of  simulation  runs. 
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Comparison  between  the  composite  curves  of  the 

two  series  of  simulations  runs 


Fig.  14  Comparison  between  the  CCs  obtained  in  the  two  series  of  simulations  with  a 
reformer  pressure  of  30  bar. 
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CONCLUSI ONS 


•  Method  for  the  design  improvement  of  a  plant  configuration  (named 
“basic")  in  which  heat  is  transferred  in  a  “black-box"  without  having  the 
heat  exchanger  structure  defined 


•  The  method  keeps  the  problem  of  the  structural  optimization  of  the 
system  independent  of  the  structure  of  the  HEN 

•  When  an  optimization  problem  is  set  up  for  the  basic  plant 
configuration,  the  most  efficient  heat  transfer  in  the  black-box  is  directly 
found  by  the  optimization  algorithm 

•  If,  conversely,  a  design  improvement  procedure  is  performed  through  a 
series  of  simulation  runs,  the  optimum  matching  between  thermal  flows 
in  the  black-box  is  found  at  each  step  of  the  procedure  using  a  sub- 
optimization  procedure  that  maximizes  the  heat  transfer  between  hot 
and  cold  flows 


On  the  synthesis  of  complex  energy  systems  with  several  heat  transfer  interactions 
by  Andrea  Lazzaretto 


CONCLUSI ONS 


•  The  method  improves  the  comprehension  of  the  heat  transfer  internal 
processes 

->  it  helps  locate  cold  thermal  flows  to  be  heated  and  properly  match 
them  with  all  the  hot  thermal  flows  available 


->  it  helps  minimize  temperature  differences  in  the  heat  transfer 
processes,  and  heat  rejected  to  the  environment  ->  maximize  plant 
efficiency 

->  it  helps  find  directions  for  design  improvement  (new  values  for  the 
design  variables  to  be  fixed  as  parameters) 
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CONCLUSI ONS 


APPLI CATI  ON  TO  THE  CROC  PLANT 

•  High  reformer  pressure  (30  bar)  with  high  methane  conversion  ratio 
(95%) 

•  Expander  to  recover  the  pressure  content  of  the  hydrogen  sent  to  the 
reformer 


•  Saturator  to  lower  NOx  production  and  improve  thermal  recovery  at  low 
temperature 


Further  advantages  by  the  water  recovery  section  because  of  the  heat 
recovered  in  the  black-box  at  low  temperature 
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Integrity  -  Service  -  Excellence 


Need  For  Energy  and  Thermal 
Management  Program  Planning 


*  Current  And  Future  Air  Force  Weapons  Systems  Lack  The 

Necessary  Power  and  Cooling  Capacity  To  Provide  Full  Systems 
Level  Capability  y  )  ...  ^  i 

*  Cooling  Capacity  Of  Fuel  Already  Fully  Utilized  Leaving  Little 
Room  For  Additional  Cooling  Needs 

*  Increasing  Speed,  Power,  And  Miniaturization  Of  Future 
Systems  Continue  To  Stress  Any  Thermal  Management 


Focused,  High  Level  Initiative  Is  Needed  To  Support 
Today’s  And  Tomorrow’s  Weapon  Systems 


Need  For  Energy  and  Thermal 
Management  Program  Planning 


*  Define  Focused  Technology  Transition  Capability 

—  Enhanced  Systems  Engineering  Process  Implementation 
—  Better  “Rules  And  Tools”  for  System  Level  Assessments 
—  Understand  Impact  Of  Technologies  on  System  Performance 

*  Apply  Systems  Engineering  Principles  --  Thermal  Management 
Considered  Early  in  System  Design 

—  Assure  Thermally  Efficient  Designs 

—  Understand  Power  and  Heat  Generation  Impacts 

*  Generate  “Game  Changing”  Approaches  To  Thermal  Management 

*  Develop  Innovative  and  More  Effective  Material  Systems 

—  Fuels,  Lubes  And  Seals 
—  Power,  Electronics,  Actuation 

Solution  Requires  Both  Better  Assessment 
Capabilities  and  More  Robust  Technologies  3 


Workshop  Goals  and  Expectations 


Address  the  Big  Questions  -  How  Much,  How  Bad 

9  What  are  the  stressing  drivers 

9  Focus  first  on  directed  energy  on  airborne  platforms 

Guide  Technology  Development 

9  From  your  perspective,  what  are  some  key  energy  and  thermal 
management  technologies  for  challenges  such  as  directed  energy 
weapons  on  airborne  platforms? 

Provide  Information  to  Guide  Resource  Allocation, 
Technology  Choices,  Technology  Designs,  and 

Research  Needs 

9  What  are  highest  payoff  areas  for  energy  and  thermal  management 
research  investments? 

9  What  systems-of-systems  analysis  methods  are  needed? 

9  How  do  we  measure  system  level  effectiveness? 
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Col  Michael  B.  Leahy,  Jr. 

Material  Group  Director 
Air  Vehicles  Directorate 
Air  Force  Research  Laboratory 


“We  give  the  Air  Force  its  Wings  ” 


Develop  and  transition  superior  air  vehicle  technology 
solutions  that  enable  dominant  military  aerospace  vehicles 


Air  Vehicles  Directorate 


Aerospace  Systems  Integration 


Control  Sciences 


Aero  Sciences 


Technology  Developer  and 
System  Integrator 


AFOSR  -  Basic  Research 


Structures 


Capability  Focused  Tech 
Investment 


Effects  needed?  When? 


S&T  needed?  When? 


What  do  we  have  to  do  to  deliver  it? 


V 


Technology  Planning  Process 

AF  CONOPS  Construct 


Air  and  Space  Expeditionary  Force 
(AEF)  Is  the  USAF  foundation 
to  prepare,  respond,  deploy, 
and  employ  for  any  task 


10  Alr  Expeditionary 


Global  Homeland 

Strike  Security  Mobility 

CONOPS  CONOPS  CONOPS 


Integration  CONOPS 
Global  Global 


Nuclear  Space  & 
Response  Response  C4ISR 

CONOPS  CONOPS  CONOPS 


^  J  Operationally  Responsive  Spacelift 

Attribute  Availability  Forecast 


BASELINE 

EELV,  Shuttle, 
Aircraft  Ops 


•  Rapid  turn  48  hrs 

•  Low  sortie  cost  $10M 

•  Vehicle  reliability  0.995 
•All  Wx  availability  90% 

•  250  Sortie  Airframe 

•  100  Sortie  Propulsion  & 
Systems 

Near  Term  (2009) 


•  Rapid  turn  24  hrs 

•  Low  sortie  cost  $5M 

•  Vehicle  reliability  0.999 

•  All  Wx  availability  95% 

•  500  Sortie  Airframe 

•  250  Sortie  Propulsion  & 
Systems 


Mid  Term  (2015) 


•  Rapid  turn  4  hrs 

•  Low  sortie  cost  $1M 

•  Vehicle  reliability  0.9998 

•  All  Wx  availability  98% 

•  1,000  Sortie  Airframe 

•  500  Sortie  Propulsion  & 
Systems 


Far  Term  (2025) 


Rapid  Turn/Low  Cost  Ops  -  High  Availability  -  High 
Confidence 

-  RLV Integration  Technologies  - 


In-house  Tech  Assessments 


A=trl 


Funding(in  $K) 

6.2  2297  2599 

6.3  20000 

SBIR  200 

Total  2297  22799 


1617  4135  3367 
30000  50000  75000 
1500 

33117  54135  78367 


Extended  Extended  Flight  Test 


50000  50000  50000  50000  50000  50000 


50300  50000  50000  50000  50000  50000 


Total  Capability  Cost:  $214,315K  based  on  $200K  class  demonstrator  +275,000K  extended  flight  tests 


Integrity  -  S e r v  i  c 


Directly  links  tech  programs  to  capability  gaps/desires 
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Long  Range  Strike 


“Ti he  capability  to  achieve  a  desired  effect(s)  rapidly 
and/or  persistently,  on  any  target,  in  any  environment, 

anywhere,  at  any  time .  ” 


•  Regional  to  Intercontinental  Reach 

•  Hold  Targets  at  Risk  for  Extended  Time 

•  Survive  &  Operate  in  Anti-Access  Environment 

•  High  Volume  and  Flexible  Payload 

•  Network  Centric  Operability 


Based  on  LRS  AoA  FY07  and  Multiple  LRS  Studies  FY01-03 


Thermal  Structures 
for  Responsive  Global  Strike 


Responsive  Global  Strike 


•  Flight  vehicle  system  heat  loads  engine, 

-  environmental  control  systems, 

-  hydraulic  systems,  and 

-  electrical  systems 

•  High  temperature  stealth  structures  - 

-  Analytical  capability  '  '  and  integrate  energy 

systems 

-  Inn*'*!  Jig  ■  •  cooled 

_i  A  system  level  approach  heir 
1  to  thermal  management 
‘  No  is  required 


Persistent  ISR 


■ 


Detect  and  track  every  target 
in  the  battlespace... 
anywhere,  at  any  time 


Based  On  ICRRA  and  GH  ORD/Annex 


dbU  degree  coverage 
On  Station  24/7L"  pS.  ""  1 
Multi-function/Multi-senSof 

Bihw  \Y\  i  jwpfri iii  i,j  ■ 


Foliage  penetration' 
E Integrated/Fused  da 
•^Penetrating  \ 


A 

Sensor  Mission  Power  Requirement 


Vol.  ft3  Wgtlbs 

Pow. 

Kw 

Radar 

0.73 

820 

62 

Radar 

0 

2000 

75 

Radar 

0 

100 

1.2 

IFF 

0 

50 

1 

ELINT 

0 

500 

0 

COMINT 

0 

0 

0 

E/O 

12 

450 

3.5 

Comms 

0.8 

52 

1 

Receiver 

Subsystem 

18 

533 

10 

Processor 

Subsystem 

7.8 

192 

4.8 

TOTAL 

39.3|. 

5000^ 

4.60-200 

For  comparison,  Global  Hawk  Integrated  Sensor  Suite 
(ISS)  weighs  about  900  lbs,  and  requires  6.5  kW  of  power 
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SensorCraft  Profile 


70,000  ft 
60,000  ft 
50,000  ft 


N.  Egress 

Cruise-Climb 

Ingress _ _ _ 

30+  hour 

3000  nm 

Loiter 

3000  nm 

>/Climb 
<200  nm 

Descent* 
<200  nm 


60,000  ft 
50,000  ft 


1  hour  Loiter 
SL 


Airframe  Around  Sensors 


Persistent  Precision  Strike 


Persistent  delivery  of  lethal  and  non-lethal 

precision  effects 


ony  Duration 


Integration) 


AADS  System 


5-10  MW  electrical 


Si 


—  4  kW  Input  Power  for  Cryogenic  Cooler  (30K)  ~1  kW 
Input  Power  for  Cathode  Heater 


<50%  efficient 


Sj 


m 


L  ) 


Idea  duty  cycle  is  100%  availability 


1-2  kHp  in 
~10%  efficient 
1000  W/cm2 


•Seconds  between  shots 
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Defending  America  by  Unleashing  the  Power  of  Innovative  Aerospace  Technology 


Energy  and  Thermal 
Management  Executive  IPT 
and  Technical  Advisory 

Group 


ETM  Workshop  Outbrief 


. .  .Enabling  Airborne 
Directed  Energy 
Weapons... 


Integrity  -  Service  -  Excellence 


Major  Themes 


*  System  Level  Assessments  (Must  Do  Now) 


—  Focus:  Get  arms  around  three  concepts  -and  hold 
true  to  them 

*  50  kw  SSL  on  a  Mid-Sized  Airborne  Platform 

*  100  kw  SSL  on  a  Tactical  Platform 

*  Airborne  Active  Denial  on  a  Mid-Sized  Airborne  Platform 

—  Bound  Problem  to  create  realm  of  possible,  not 
preselect  point  solution 


Sim  Based  R&D 

The  Engine  in  Capability  Planning 


Conceptual 

Studies 


•  Vehicle  Concepts 

•  Technology  Assessment 

Sensitivities  I  pPg 
studies  I®:. : 


7 

£ 

/ 

/ 

i 

Tech  Programs 


Technology 
High  Fidelity  Trades 
Analysis 
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Current  Status 


Thermal  and  Energy  Management  Being 
Worked  At  Component  And  Subsystem  Level 


Shortfall 

*  Methods  And 

Capabilities  To  Rapidly 
Understand  The  Impacts 
Of  Technologies  At  The 


System/Capability  Level 


Assessment 


Identified  Key  Technology 
Development  Drivers 
What  is  Required  to. . . 

•  Fly  how  Long 

•  Shoot  so  Far 

•  Carry  How  Much 

n»  Power  Generation 
•  Device 

•  Thermal  Tech  Trades 
•  Power  and  Thermal  Storage 


Mission  Definition  and 
Analysis 

Effectiveness  and  LCC 
Assessments 

Assessment  of  Technology 
Impa"*" 


Linkage  of  KPPs  to 

Capability 

Link  to  Technology 

Demonstrations 


High  Fidelity  Cooling  Analysis 
Structural  Analysis 
Sophisticated  Device  Modeling 
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Major  Themes  (Con’t) 


*  Theme:  Increasing  Weapon  Efficiency  will  only  help 
the  thermal  management  issues,  not  solve  the  issue 

*  Theme:  Need  for  advanced  fluids  and  materials 

*  Theme:  Mature  Two-Phase  Flow 

*  Theme:  Energy  Harvesting  (Way  Cool  If...) 

*  Theme:  Energy  Storage  Systems 

*  Phase  Change  Materials  (PCM),  Adsorption, 
Chemical  reaction,  ... 


*** 


Energy  and  Thermal  Management 
Executive  IPT  Vision 


Includes  all 
players  & 
schedules 


The  systems 

approach 

_  J 


The  pipeline 

The  legacy 

_ _ _ / 


To  develop/promote  &  establish  a  robust 
energy  and  thermal  management  research  & 
development  program  that  is: 

*  Integrated  across  AFRL  -  Maximizes  use  of 
other  organizations’  efforts 

*  Focused  on  capabilities  at  the  system  level 

*  Provides  for  a  balance  of  6.1 ,  6.2  and  6.3 
research  &  development 

...to  unleash  the  power  of  innovative 
technology. . . 


Cross  Cutting  Executive  IPT  Formed  and  Working 

Integrated  Approach 


Air  Force  Research  Laboratory 


Defending  America  by  Unleashing  the  Power  of  Innovative  Aerospace  Technology 


Energy  ancPThermal 
Management  Stressing 
Scenarios 

An  Air  Vehicles  Directorate  Point  of  View 


David  M  Pratt,  PhD 


OR 


Integrity  -  Service  -  Excellence 


/.ir!;r,rr.r 
Dental  Syrir, 


Tactical  Laser 

Mid-Sized  Airborne  Platform 


MtMrfwcc  J1 


TMS  Requirements 


•1-2  kHpin 
•~10%  efficient 
•1000  W/cm2 


Seconds  between  shots 


